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Abstract: Holographic models in the T = universality class of QCD in the limit of large 
number N c of colors and Nf massless fermion flavors, but constant ratio xt = Nf/N c , are 
analyzed at finite temperature. The models contain a 5-dimensional metric and two scalars, 
a dilaton sourcing TrF 2 and a tachyon dual to qq. The phase structure on the T, xj plane 
is computed and various 1st order, 2nd order transitions and crossovers with their chiral 
symmetry properties are identified. For each x/, the temperature dependence of p/T 
and the condensate (qq) is computed. In the simplest case, we find that for Xf up to the 
critical x c ~ 4 there is a 1st order transition on which chiral symmetry is broken and the 
energy density jumps. In the conformal window x c < Xf < 11/2, there is only a continuous 
crossover between two conformal phases. When approaching x c from below, Xf — > x c , 
temperature scales approach zero as specified by Miransky scaling. 
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1. Introduction 

QCD in the Veneziano limit [1], 

Nf 

N c — > oo, Nf — > oo, T7" =x / fixed, A = gy M N c fixed, (1-1) 

is expected to display a host of interesting and mostly non-perturbative phenomena, in- 
cluding: 

• The "conformal window" with a nontrivial infrared (IR) fixed point, which extends 
from Xf = to smaller values of xj. The region xj — > ^ has an IR fixed point 
while the theory is still weakly coupled, as was analyzed in [2] (see also [3]). 

• It is expected that at a critical Xf = x c , the conformal window will end, and for 
Xf < x c , the theory will exhibit chiral symmetry breaking in the IR. This behavior 
is expected to persist down to Xf = 0. For Xf>x c the IR theory is a conformal field 
theory at strong coupling, that progressively becomes weak as Xf — > 4r. 

• Near and below x c , there is the transition region to conventional QCD IR behavior. 
In this region the theory is expected to be "walking": This means that the theory 
appears to be approaching the IR fixed point as the coupling evolves very slowly for 
many e-foldings of energies. But chiral symmetry breaking is nevertheless triggered 
and in the deep infrared the coupling diverges as in QCD. The slow evolution of the 
coupling has been correlated with a nontrivial dimension for the quark mass operator 
near two, rather than three (the free field value). IR observables are expected to obey 
the Miransky scaling [4, 5] as Xf — > x c from below. 
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• New phenomena are expected to appear at finite density driven by strong coupling 
and the presence of quarks. These include color superconductivity [6, 7]. In this case, 
however, gauge invariant vevs are effectively double trace operators and the phase 
structure is determined at the next to leading order in 1/N C . 

The existence of the "walking" region makes the theory extremely interesting for appli- 
cations in dynamical electroweak symmetry breaking (technicolor). This has also motivated 
an intensive lattice Monte Carlo work during recent years [8, 9, 10]. The bulk of this work 
has been done at zero temperature; recently there appeared the first attempts to go to 
finite T for QCD with N c = 3, N f up to 8 [11, 12, 10] and for non-QCD-like theories [13]. 
Chiral effective theories have also been applied [14, 15, 16, 17, 18, 19, 20, 21]. 

The aim of the present work is to study a class of holographic bottom-up models (V- 
QCD) that belong to the universality class of QCD with massless quarks in the Veneziano 
limit [22] at finite temperature and zero chemical potential. We will calculate the tem- 
perature dependence of the free energy density (= —pressure = —p(T)) and of the quark 
condensate ({qq)(T)). The former acts as an effective order parameter for deconfinement 
(at large N c ), for which there is no true order parameter associated with a symmetry. 1 
The quark condensate is a true order parameter for chiral symmetry if quarks are massless. 
The calculation is carried out for the full range of x/, < xt < 11/2. 

Discontinuities or rapid variations in pressure (or energy density) and quark condensate 
can be used to define phase boundaries associated with deconfinement and chiral symmetry 
restoration temperatures T^xj) and T x (xf). We will use the usual nomenclature: If the 
nth derivative of p(T) is discontinuous, the transition is of nth order. We also consider 
continuous crossovers which are identified by using the scaled quantity (e — 3p)/T 4 . Its 
maximum defines the crossover temperature T crossover (a;y). The phase diagram is defined 
as a plot of all phase boundaries on the (xf,T) plane. The phase diagrams we present will 
also contain a rich structure of metastable states, namely local (but globally subleading) 
minima of the free energy. 

In the holographic approach the thermal transitions will be transitions between various 
5-dimensional black hole and "thermal gas" metrics and the nomenclature of transitions, 
explained later in great detail, will be correspondingly different. The holographic approach 
is constrained but not fully constrained and we cannot give a precise prediction of the 
phase diagram of hot V-QCD. We can state the most plausible behavior but we can also 
mention a few other alternatives. We will always find the analogues of and T x , but we 
will also find transitions with no obvious QCD interpretation. Whether these reflect real 
physics of hot QCD in the Veneziano limit or whether they are artifacts of the holographic 
approach will be an interesting problem for further study. 

The usual expectation is that there is a 1st order line at = T x ; in the large iV c 
limit one can actually prove that T x > [14, 15]. The main class of our predictions 
reflect these properties: for smaller Xf we find that deconfinement and chiral symmetry 
restoration coincide, but for Xf approaching x c the deconfining and chiral transitions can 
become separate so that T x > (see, for example, Fig. 13 below). The chiral transition 

1 A related one, used commonly in lattice work, is the expectation value of the Polyakov loop. 
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is then of 2nd order (and mean field type). Furthermore, for smaller Xf the separate 2nd 
order chiral transition is in the metastable region so that it can be reached if the system 
is supercooled [23]. One might here add that T x < for stable phases may be reached at 
large chemical potential [24, 25]. 

The starting point of our finite temperature analysis is the T = holographic model 
introduced in [22], based on previous theoretical ideas in [26, 27, 28, 29, 30]. Moving to 
finite T implies studying black hole solutions of the action in [22] . A defining characteristic 
of this class of models is that they contain full backreaction between the duals of the color 
and flavor degrees of freedom. Earlier work [31, 32, 33, 34] on thermodynamics in such 
bottom-up models imposed quasiconformality directly on the beta function of the theory. 
One should note that walking behavior and the related "conformal transition" at xj = x c 
have also been studied in top-down models [35, 36, 37, 38], as well as in simpler bottom-up 
models [39, 40] which do not attempt to model the backreaction. See also the review [41] 
on introducing backreacted flavor in the top-down models. 

In this introduction we will first describe the special properties of V-QCD from [22] 
and then discuss general properties of its black hole solutions. Section 2 will contain a 
detailed discussion of the action of the model and of the two characteristic classes of scalar 
potentials. Section 3 presents the Einstein equations of the model, describes how they are 
numerically solved and, finally, how thermodynamics is computed from the numerical bulk 
fields. A particularly delicate issue here is the fixing of the quark mass m to zero. We also 
briefly comment on fixing m to some nonzero value. An extensive list of numerical results 
is given in Section 4. From these, the types of phase transition lines the models predict 
are determined. In Section 5, techniques for computing the condensate are described and 
several numerical results are given. One should note that this, as well as many other 
numerical issues in the model, are technically very demanding. Finally, Section 6 contains 
a discussion of what are the effects of making the quark mass nonzero. Several detailed 
considerations are collected in Appendices. 

1.1 V-QCD at zero temperature 

In [22] a class of bottom-up holographic models was introduced (named V-QCD) and 
shown to be in the universality class of QCD in the Veneziano limit at zero temperature 
and density. These were 5-dimensional models of two scalars coupled to gravity. One of 
the scalars, the "dilaton" A, is dual to Tr[F 2 ] (the QCD gauge coupling constant, or more 
precisely the 't Hooft coupling). The other scalar, the "tachyon" r, is dual to the quark 
mass operator qq. The potentials and interactions were modeled along successful bottom- 
up models for YM, namely Improved Holographic QCD (IHQCD) [26, 27, 28] and the idea 
that string theory tachyon condensation describes chiral symmetry breaking [29, 30, 42]. 

The bulk action considered was 



S = Sg + Sf 




(1.2) 
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with A the 't Hooft coupling (exponential of the dilaton <fi) and the tachyon 2 action 3 

S f = -x f M 3 N* J d 5 x V f (X, t) ^ det( 9flu + K^df.rduT). (1.3) 

The pure glue potential V g has been determined from previous studies [27]. The 
tachyon potential Vf(X, r) must satisfy some basic properties determined by the dual theory 
or by general properties of tachyons in string theory: (a) To provide the proper dimension 
for the dual operator near the boundary (b) To exponentially vanish for r — > oo. The 
function k(X) captures, among other things, the transformation from the string frame to 
the Einstein frame in five dimensions. The class of potentials that were investigated in [22] 
are of the form 

y / (A,r) = Vo/(A)e^ A ^ 2 . (1.4) 

In the Veneziano limit, the back-reaction of the flavor sector on the glue sector is fully 
included. 

As with IHQCD, it was arranged that the theory is asymptotically AdS in the UV up 
to logarithmic corrections in the bulk coordinate. The function Vo/(A) is such that the 
potential V g (X) — XfVof(X), when the tachyon has not condensed (r = 0) has an extremum 
4 at a finite value A = A*. As we approach the Banks-Zaks region [2], xf — > 4r, the value of 
A* approaches zero. Without the tachyon, r = 0, the equations of motion imply that also 
/3(A*) = 0, i.e., A* is an IR fixed point. When the dynamics of r is included, the system 
approaches A* but is driven away from it as long as xj < x c (see Fig. 7 of [22]). 

The dimension of the chiral condensate was calculated in the IR fixed point theory from 
the bulk equations. It was found that it decreases monotonically with Xf for reasonably 
chosen potentials. It crossed the value 2 at xt = x c where x c corresponds to the end of the 
conformal window as argued in [43]. 

The lower edge of the conformal window x c lies in the vicinity of 4. Requiring the 
holographic /3-functions to match with QCD in the UV, we find that 

3.7<x c <4.2, (1.5) 

where the bounds are not strict but hold approximately for potentials that have smooth 
A-dependence in the UV. 

There is also a phenomenological heuristic argument for the value x c ~ 4, simply from 
counting degrees of freedom. At low T chiral symmetry is broken and the massless degrees 
of freedom are Nj Goldstone bosons. At large T there are 2N% + | N c Nf weakly coupled 
degrees of freedom. These numbers are equal for Xf = 4. Conformal window and the 
location of its edge was also discussed within holographic frameworks related to V-QCD in 
[44, 33]. 

2 We have taken the tachyon to be real and diagonal in flavor space. 

3 To find the vacuum (saddle point) solution we have set the gauge fields A^' R dual to the QCD currents 
to zero, as they are not expected to have vacuum expectation values at zero density. We have also suppressed 
the Wess-Zumino terms as they also do not contribute to the vacuum solution. 

4 The extremum may exist for all < Xj < ^ or may disappear at some small Xf. No changes in the 
phase structure at zero temperature for these two cases were found in [22]. 
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Apart from x f , there is a single parameter in the theory, namely Aqq D wnere 171 is the 
UV value of the (flavor independent) quark mass. For each value of Xf, the bulk equations 
were solved with fixed sources corresponding to fixed m, Aqcd- The vevs were determined 
such that the solution is "regular" in the IR. The notion of regularity is tricky even in the 
case of IHQCD (pure glue), as there is a naked singularity in the far IR. For the dilaton 
this has been settled in [26, 27]. For the tachyon the notion of regularity is different and 
has been studied in detail in [30]. 

The regularity condition was implemented in the IR. After solving the equations from 
the IR to the UV (this was done mostly numerically), there is a single parameter that 
determines the solutions as well as the UV coupling constants and vevs, and this is a real 
number To controlling the value of the tachyon in the IR. This number reflects the single 
dimensionless parameter — of the theory. 

For different values of Xf and m the following qualitatively different regions were found 
in [22]: 

• When x c < Xf < 11/2 and m = 0, the theory flows to an IR fixed point. The IR 
conformal field theory is weakly coupled near Xf = Q and strongly coupled in the 
vicinity of x c . Chiral symmetry is unbroken in this regime (this is known as the 
conformal window). 

• When x c < Xf < 11/2 and m ^ 0, the tachyon has a nontrivial profile, and there is 
a single solution with the given source, which is "regular" in the IR. The IR theory 
is a theory with a mass gap. 

• When < xj < x c and m = 0, there is an infinite number of regular solutions with 
nontrivial tachyon profile, and a special solution with an identically vanishing tachyon 
and a nontrivial IR fixed point. The infinite number of solutions with nontrivial 
tachyon are classified by their number of zeros. The solution with the lowest free 
energy is the one with no zeros. 

• When < Xf < x c and m ^ 0, the theory has vacua with nontrivial profile for the 
tachyon. For every non-zero m, there is a finite number of regular solutions that 
grows as m approaches zero. 

In [22] two large classes of tachyon potentials were identified. Potentials in class I, 
have a(A) constant in (1.4). In this case the tachyon diverges exponentially in the IR for 
the regular solution 

t ~ roexpfCr], (1-6) 

r— >-oo 

where C is a known constant (see Appendix B) and tq is the only integration constant 
controlling the solution. It determines the source (mass) in the UV. Potentials in class II, 

2 

have a(A) ~ As as A — > oo, and a tachyon that diverges in a milder way in the IR as 



rsJ 



C^r^n, (1.7) 
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where again C is known and r\ is the single integration constant controlling the regu- 
lar solution. The qualitative conclusions above and below were valid for both classes of 
potentials. 

In the region Xf < x c where several solutions exist, there is a interesting relation 
between the IR value tq controlling the regular solutions, and the UV parameters, namely 
m. This is determined numerically, and a relevant plot describing the relation between m 
and ro at fixed Xf is in Fig. 1. 

The solutions are characterized by the number of times n the tachyon field changes 
sign as it evolves from the UV to the IR. For all values of m there is a single solution with 
no tachyon zeroes. In addition, for each positive n there are two solutions 5 which exist 
within a finite range < m < m n , where the limiting value m n decreases with increasing 
n, and one solution for m = 0. In particular, for large enough fixed m, we find that only 
the solution without tachyon zeroes exists. 

For m 7^ 0, out of all regular solutions, the "first" one without tachyon zeroes has 
the smallest free energy. The same is true for m = 0, namely the solution with nontrivial 
tachyon without zeroes is energetically favored over the solutions with positive n as well 
as over the special solution with identically vanishing tachyon, which appears only for 
m = and would leave chiral symmetry unbroken. Therefore, chiral symmetry is broken 
for Xf < x c . 

In the region just below x c , [22] found Miransky scaling for the chiral condensate. As 

Xf r X c , 



y y x c Xf I 

For x > x c , let mi^(x) be the mass of the tachyon at the IR fixed point and £ir(x) the IR 
AdS radius. The coefficient K is then fixed as 



The behavior at and below the conformal transition at Xf = x c is to a large extent 
independent of the details of the model. In particular, no information on the nonlinear 
terms in the tachyon EoM is needed or how the IR boundary conditions are fixed. In the 
same region, "walking" of gauge coupling is realized. The YM coupling flows from small 
values to values very near A*, remains approximately constant for many e-foldings of energy 
(in this regime the tachyon remains small), and then runs off to infinity, driven by a large 
value of the tachyon field in the IR. The walking is related to a long section of the solution 
which is similar to the one studied in earlier bottom-up models for walking [39]. 

The finite temperature analysis of V-QCD amounts to studying all black hole solu- 
tions with appropriate boundary conditions. To start with, any zero temperature solution 
becomes a candidate saddle point at finite temperature by compactifying time on a cir- 
cle of radius f3. Any other competing black hole solution must have the same boundary 
conditions as well as a regular horizon in the IR. 

As m and —m are related by a chiral rotation by 7r, we can take m > 0. 




(1.8) 



K = 



7T 



(1.9) 
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As the dilaton always has a nontrivial UV source, it will always have a nontrivial 
profile in the black-hole solutions. With the tachyon, things can be different. In the 
massless case, its source is zero. Therefore there are two possible options (as in the zero 
temperature configurations discussed above): either it is identically zero (if the vev (qq) is 
also zero) or it is non-zero (implying a non-zero vev). 

Therefore we have two large classes of black holes in the massless case: those with 
r = and those with We will first consider the tachyon-free class. 

1.2 Black holes without tachyon hair 

If r = 0, we have black holes in a single scalar theory, with potential V(A) = V g (\) — 
XfVof(X) from (1.4). This is a potential with an extremum for xt ^ 6 and no extremum 
when Xf = 0. 

Black hole solutions for such potentials were discussed in generality in [27]. After fixing 
all invariances, they are characterised by a single IR constant, A^, the value of the dilaton 
at the horizon. The plot relating the temperature T to A^ contains important information 
about the thermodynamics of such black holes. Small values of Xh denote large black-holes 
whereas larger values of A^ correspond to smaller black holes (smaller horizon size and 
entropy). In all plots of this paper, dilatonic black holes without tachyon hair are denoted 
by red lines in the respective (A^, T)-diagrams, and we shall call the corresponding function 
T u (X h ). 

When Xf = 0, A can become arbitrarily large at zero temperature, implying that A/j 
can also be arbitrarily large for the finite temperature configurations. At finite temperature 
there are two branches: large black holes which are stable and small black holes which are 
unstable. If T^(A^) < the black-hole thermodynamics is stable, otherwise it is unstable. 
There is a minimum temperature above which such black holes exist as shown, for example, 
by the black line in Fig. 22 (left or right). 

When Xf > 0, we have two possibilities. The first is that the potential V e s(X) = 
Vg(\)—XfVof(X) has an extremum at A — > A*(x) for all < xj < ^, with A* (xf W) — > 
and \*(xf — > 0) — > oo. The second is that such extremum only exists for x > x*, where 
x* < x c . We shall denote these potentials with a star subscript. 

At finite temperature, and when the potential V e s has no extremum, the black hole 
without the tachyon hair exists for all positive A^. For the potentials studied here, function 
T u (Xh) is qualitatively similar to the YM case (xf — > 0) [27]. As shown in Fig. 17 (top-left) 
and in Fig. 19 (left), there are two black hole branches, which exist above some minimum 
temperature. The branch at low A^ is thermodynamically stable, while the large-A/j branch 
is unstable. 

When the extremum is present, < A^ < A* (a;). The temperature of the black- hole 
corresponding to Xh = A*(x) is T = 0, while that of Xh — > has T — > oo. There is no 
minimum temperature here. For any temperature there is always at least one black-hole 
solution. There are several possibilities that are shown as red lines in Figs. 7 (left), 9 (top), 
10 (left) and 12 (left). 

6 The extremum may also exist only for Xf above some fixed x*, see the discussion further below. 
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When Xf is large, but still smaller than x c , the T = T u (Xh) relation is one-to-one but 
contains a bump (a change of concavity) like in Fig. 9 (top). Then this is accompanied by 
a crossover behavior, signaled by a bump in the trace of the stress tensor (e — 3p)/T 4 , (aka 
interaction measure) as shown in Fig. 9 (bottom-right). 

At low enough Xf, the relation T = T u (Xh) is not always one-to-one, as can be seen in 
Fig. 10 (left) or in Fig. 22. Then there are points where T^A/J = 0. In such a case there 
can be a first order transition between the stable branches of the black hole solutions. This 
is a remnant of the deconfining transition at Xf = (pure YM). In Fig. 22 both left and 
right several curves in the (T, A/J-plane for different xj indicate the successive structure of 
dilaton black holes (red lines). The black line corresponds to the pure YM (xf = 0) limit. 

When x > x c we are in the conformal window. The only black holes that exist here are 
those without tachyon hair. The relation T = T u (Xh) is monotonic and there is a continuous 
transition to the black- hole phase at T = + , as in the AdS case in the Poincare patch. The 
thermodynamic functions, especially the interaction measure, show a crossover maximum 
at a temperature that is moving towards the UV as Xf — > . 

1.3 Black holes with tachyon hair and zero quark mass 

When r / we have black holes in the two scalar theory. The tachyon starts as ~ r 3 
near the UV boundary as the source (quark mass) vanishes. In all plots of this paper, 
such black holes (with both dilaton and tachyon hair) are denoted by blue lines in the 
respective (A^, T)-diagrams, and we shall denote the corresponding functions by Tb(A^). 
They are still one parameter solutions and can be parametrized again by the value A^ of A 
at the horizon, which translates into the temperature. These black holes usually exist for 
all Xf g]0,x c [ and our discussion below focuses in this region. 

Because the presence of the nontrivial tachyon perturbs and annuls the possible non- 
trivial IR fixed point, for such black-holes, A/j can take arbitrarily large values. This can 
be seen from the blue lines in Figs. 7 (left), 9 (top), 10 (left) and 12 (left). For all such 
black holes, the chiral condensate is determined by the regularity of the black hole solution. 
It decreases as A^ decreases, and at some point it vanishes. At this point, the blue line 
in the (A^, T)-diagram merges with the red line corresponding to a A^ that we call A cnc j 
throughout the paper. This can be seen in all the figures mentioned above. 

The shape of the blue line can vary as a function of Xf and the type of potential. There 
are three typical examples of shapes: 

• Simple lines that are monotonic as the one depicted in Fig. 12 (left). This is an 
example of a monotonic blue branch where all such black-holes are thermodynamically 
unstable. Moreover, they have a minimum temperature. In such a case, they can 
never be thermodynamically dominant. At some temperature the vacuum thermal 
solution is dominated by a dilaton black hole on the red line, and the chiral restoration 
transition is 1st order. 

• Lines with two branches as the one depicted in Fig. 10 (left). Here the blue line has 
two parts one (to the left) that is thermodynamically stable and another to the right 
that is thermodynamically unstable. In such a case, the system is in the thermal 
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m 

^"^^^^ T ° 

Figure 1: Plot of the UV Mass parameter m, as a function of the IR scale To in (1.6), for xj < x c . 
The vertical solid blue line marks the end-point of the existence of regular solutions. The dashed 
red line indicates the appearance of more than one regular solution with the same value of m. 

vacuum solution at low enough temperatures, then jumps with a 1st order transition 
to the tachyon-hairy solution (the part of the blue line that is thick in Fig. 10 (left)) 
which still break chiral symmetry, and then eventually smoothly transits to the red 
line at the point where the blue and red lines merge, via a chirally-restoring 2nd order 
transition. 7 

• Lines with more than two branches as the one depicted in Fig. 11 (left). In this 
example the blue line has four branches, two unstable and two stable. There are 
in total three phase transitions here, first from the vacuum thermal solution to the 
rightmost blue thick branch, then to the intermediate thick blue branch and finally a 
2nd order chirally restoring transition to the red branch at the point they touch. In 
this case there are two 1st order transitions between three chirally breaking phases, 
and a 2nd order one to the chirally symmetric phase. 

A concrete overall view of the Xf dependence is presented in Fig. 2, in which T(Xh) 
is plotted for potentials of type II with SB normalisation (definitions specified later) for 
various xt. One sees clearly how the pure (black) YM curve is approached for xt — > 0. 
The thick curves represent stable phases; when a thick curve ends, the system makes a 1st 
order transition to the low T phase. When thick curves change from red to blue curves, a 
2nd order transition to a chirally broken phase takes place. For a more accurate picture of 
small Xf, see Fig. 22. 

1.4 The phase structure of different V-QCD models 

There are three main ingredients that characterize a priori different QCD models which, 
however, have the same phase structure and qualitative behavior at zero temperature: 

• The asymptotics of the tachyon solution in the IR. This is controlled by the behavior 
of the function a(A) in the tachyon potential in (1.4). When a(A) is constant, the 

r It may also happen that the thermodynamically stable branch is only metastable, in which case the 
system jumps directly to the black hole branch without tachyon hair, and chiral symmetry is thus restored 
at this 1st order transition. The more complicated branch structure discussed in the next bullet may 
similarly contain metastable branches. 
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Figure 2: T(A/ l )/A plotted for potentials of type II with SB normalisation (definitions specified 
in text) for various xt marked on figure. Thick curves represent stable phases. For more details at 
small Xf, see Fig. 22. 

tachyon diverges exponentially in the IR, and we call such potentials of type I. When 

2 

a(A) diverges as As in the IR (A large) then the tachyon diverges as a square root in 
the IR, and we call such potentials of type II. 

• For any potential, the UV constant factor Wo of Vo/(A) in (1.4), defined in (2.13) can 
vary in finite range, which in appropriate units is ]0, as in (2.21). We pick for each 
type of potential three indicative values of CqWq that in general might give different 
physics, namely 0, y|, and j|. 8 We also consider xj-dependent value, specified in 
(2.23) that corresponds to the normalization of the UV degrees of freedom of the free 
energy to the Stefan-Boltzmann limit of QCD. 

• A final variation can be obtained on all of the above by using a glue potential V^(A) — 
XfVof(X) in (1.4) that has 

(a) an extremum for all xj in the appropriate range, xj €]0, 

(b) an extremum only in part of this range, x* < xj < We will denote the 
potentials in this case by a star subscript. 

According to the above options PotI*(W / o = 0) denotes a potential in the type I class, 
with Wq = and an IR critical point that exists only down to a finite X*. 

Let us then summarize the phase structure of the model as Xf and the temperature 
are varied (at zero quark mass). In general one expects the phase diagram of Fig. 8, so 
that for < Xf < x c there is the 1st order transition at finite temperature, which also 
separates the chirally symmetric and broken phases. For Xf > x c the low temperature and 

8 Notice that the exactly zero value of Wo is actually excluded, because it predicts wrong anomalous 
dimensions for quark mass or the chiral condensate [22]. We anyhow consider it as the limiting case of the 
allowed solutions. Moreover, Wo may exceed the upper limit of 24/11, if Xf dependence is allowed. 
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high temperature configurations correspond to a tachyonless black holes, and, one expects 
a continuous crossover between these two. 

For the various potentials presented above, this phase diagram is indeed obtained in 
the zeroth approximation, but for Xf < x c there are additional details which depend on 
the choice of potentials as follows. 

• For potentials I the phase structure depends strongly on the choice for Wq (see 
Fig. 18). For the lowest value Wq = 0, there is only one 1st order transition at 9 
T = Td for all < Xf < x c , except possibly for Xf very close to x c , where solving the 
phase diagram numerically becomes demanding. As Wq is increased, a complicated 
structure appears near Xf = x c , where we have two 1st order transitions between chi- 
rally broken phases, and the restoration of chiral symmetry at a 2nd order transition 
at even higher temperature. At even higher Wq the 1st order transitions combine 
again into a single one, but the separate 2nd order transition continues to exist for Xf 
close to x c . At low Xf, there is also a surprising change as Wq increases. The chiral 
symmetry breaking phases disappear, but there is a 1st order transition between two 
chirally symmetric black hole phases at a finite temperature instead. 

• For potentials II the dependence on Wq is milder (see Figs. 13 - 16). At high Wq, for 
low Xf up to some value x x , there is only the 1st order transition at 10 T = T^. When 
x x < Xf < x c , the chiral symmetry restoration takes again place at a 2nd order 
transition at T x such that T x > Tj. For decreasing Wq, x x increases, and finally 
disappears by joining with x c . 

• For the potentials I*, the phase structure is the standard one for high Xf, i.e., a 2nd 
order transition and a 1st order one with critical temperatures T x > within a range 
x x < Xf < x c , with the former separating the chirally symmetric and broken phases 
(see Fig. 19). For lower Xf there is only one 1st order transition. For Xf < 2, in the 
region where the effective potential does not admit an extremum, chiral symmetry 
is intact at all temperatures. We find a single 1st order transition between chirally 
symmetric thermal gas and black hole phases. 

• For potentials II*, the phase structure is simple (see Fig. 17): there is a single 1st 
order transition for all Xf s]0,x c [. In particular, the system is in a chirally broken 
phase at low temperatures, even in the region of low Xf where the effective potential 
does not have an extremum. 



2. Defining V-QCD 



2.1 Gravity action of the model 

The action of V-QCD is [22] 





°T d = T h in Fig. 18. 

°T X — r en( j of Figs. 13 - 16 when it is in the stable brach. 
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where 



11 



-9 R 



4 (d\f 



3 A 2 



Vg(X)) - V f (\, t) Vdet (g ab + k(A, r)d a r d b r) 



R+l-lg^d^d^ + VgiX)] - Vf(X,r)y/l + g rr k(X(t))t' (r) 



The metric Ansatz is 



ds 2 = b 2 (r) 



-f{r)dt 2 + dx 2 + 



dr 2 

W) 



, b(r) = e A{r) 



C 



uv 



r— >0 r 



and the two scalar functions, 1/A sourcing F 2 and r sourcing (qq), are 



A = A(r) 



r = r(r). 



(2.2) 



(2.3) 



(2.4) 



In the second form \/—g has been factored out of the DBI action. The Gibbons-Hawking 
counter term would be 



S G h = ~ J d 4 x 
with, for a hypersurface r =const, 

K = 



"7 



7(J{r) f'(r) 



(2.5) 



(2.6) 



Notice also that we have set the gauge fields ^4l,_r, which are dual to the left and right 
handed fermion currents, to zero, and neglected the Wess-Zumino terms. These terms do 
not affect the thermodynamics of the models. 

The background solution of the dilaton A(r) and the warp factor A(r) are identified 
as the 't Hooft coupling and the logarithm of the energy scale of the dual field theory, 
respectively [26]. As a matter of convention, we shall fix the normalisation of A(r) so that 
its relation to the perturbative QCD coupling is 



A(r) 



9 2 (b(r)) 
8tt 2 ' 



(2.7) 



The results of the model are independent of this normalisation, changing A — > X/Xq one 
simply has to change the potentials by V(A) —> V(AoA). The convention of [22], for 
example, is obtained by shifting by Ao = l/(87r 2 ). 

Important ingredients of the model are the relation of the bulk fields at r to the QCD 
beta and quark mass anomalous dimension functions evaluated for a coupling at scale b(r). 
Motivated by the connection to field theory, one defines 







dX 
db/b 



X\A) 



-b X 2 - M 3 - 6 2 A 4 



1 = r'(A). 



(2.8 



11 Notice that for notational simplicity we have absorbed a factor of Xf, which is visible in Eq. (1.3), into 
V f (X,T). See also Eq. (2.10) below. 
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Matching with the perturbative expansion of the QCD beta function gives 

6 = I (11- 2x f ), h = i (34 -13s/). (2.9) 

The action contains the gluonic potential V g (X) and the fermionic potential Vf(X, r), 
which will be specified to the form 

V f (X,r) = x f V f0 (X)e- a ^ T2 . (2.10) 

The detailed form of these and the functions k(A), a(X) will be discussed in the following 
subsections. 

2.2 Construction of the potentials 

The potentials can be constructed in stages. First one fixes the potentials V g (X) and 
Vfo(X) up to order A 2 in the UV, using the two scheme independent coefficients of the 
beta function. This analysis is simplified by the fact that the tachyon decouples in the 
UV. Next one fixes the UV behavior of the functions a(A) and k(A), which parametrize 
the tachyon dependence of the action using the similarly scheme independent UV running 
properties of the quark mass and the condensate. Finally, one fixes the large A behavior 
of the potentials by requiring that the model reproduces known features of QCD in the 
IR, such as confinement, linear Regge trajectories, and reasonable zero-temperature phase 
structure. We shall discuss the various steps in detail below (see also [22]). 

2.2.1 The potentials from the beta function in the UV 

In the UV, since the tachyon vanishes much faster than the dilaton, we can first set it to zero. 
Then the dilaton profile can be linked to the effective potential V e g(X) = V g (X) — XfVfo(X) 
[22] by using Einstein's equations [26]. Defining (3 = dX/dlnb = — 6qA 2 — 61 A 3 , to order A 2 , 



12 

V 9 ~ x f V of = exp 



A 



.18 



12 

T 2- 



l + §M + (|jg + ^)A a 



(2.12) 



= V - x f W + (Vi - x f Wi)X + (V 2 - x f W 2 )\, (2.13) 

where we expanded 

V g = V + V 1 X + V 2 X 2 + 0(X 3 ), V f0 = W + W 1 X + W 2 X 2 + 0(X 3 ), (2.14) 
and where we have introduced an Xf dependent AdS radius 

£ uv = C(x f ). (2.15) 
Applying equation (2.12) for Xf = we have for the gluonic potential 

y ^^ 1 + B 6 ™ A+ 23 M M)^ M M^ (216) 

12 / 88 x 4619 , 2 \ . w , 
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where bJ M are the values of b{ for xj = and Cq = C(xj = 0). In practice, one usually 
sets the (dimensionful) quantity Co = 1. 

By using equations (2.12) and (2.13) one can now solve for the coefficients of the 
fermionic potential: 

x f C 2 W = J^J, (2.18) 



_ 32 I b YM 



c 2 

00 



C 2 



12 • 8 
27 



ll-(ll-2x/)-^ n 
^uv 



(2.19) 



12 
81 
12 
729 



23( 







+ 366 1 YM -(236^ + 366i ; 

^uv 
C 2 1 

4619 - (4619 - 1714s/ + 92xJ)-^ 



uv 



(2.20) 



These equations still involve one free parameter, which can be taken to be either Wq or 
Ajv- We shall study two ways to fix this parameter. First, we can take Wo to be constant. 
In this case [22] 



< £§W < 
and the xj-dependent AdS radius is given by 



24 

TT' 



(2.21) 



£uv 



l-±C 2 W -x f 



(2.22) 



Second, we can make a special x/-dependent choice, which (as we shall show later) auto- 
matically normalises the free energy at large T to Stefan-Boltzmann: 

C vy = C (l + lx f ) 1 '\ (2.23) 

Further, we have to fix the A dependence of the functions a(A) and k(A) in the tachyon 



part 



x f V f0 (\)e~ a ^ r2 Vl + 5 rr ^(A(r))f 2 , 



(2.24) 



of the action, where g rr = f/b 2 . The leading logarithmic term to the UV expansion of the 
tachyon should be (remember that the energy dimension of r is — 1) 



__3_ 

r(r) / C\jy = mr (—In Ar) b o = mr (— In Ar) 2b o 



(2.25) 



to satisfy the scheme independent UV running of the quark mass. Here 70 = 3/2 is the 
leading coefficient of the anomalous dimension of the quark mass in QCD, 7(A) = 7oA + - • • . 
By using the tachyon equation of motion one sees that this requires that for small A, 



«(A) 
a{\) 



C 



uv 



1 



60 + I A + A + 



(2.26) 
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2.2.2 Large A behavior of the potentials 

To specify the full potential V g (X) — x fVfo(\)e~ a ^ T2 we have to continue the small A 
expansions to large A. The guideline is quark confinement and chiral symmetry breaking 
at small Xf and the appearance of an infrared fixed point at some Xf = x c (see [22]). Since 
there is no unique path to the result, we present the final forms of the potentials we use 
and motivate them. 

We use the gluonic potential 



Vg(X) 



12 

T 2 



88A 4619A 2 y/1 + ln(l + A) 
+ ~27~ + 729 (1 + A) 2 / 3 



(2.27) 



which is constructed from the expansion (2.16) by simply multiplying the A 2 term by the 
confinement factor 

Vl + ln(l + A) 

(1 + A)V3 • ( 2 - 28 ) 

Then V g has the proper large- A behavior [26] but the small- A behavior is left intact. One 
could add scale factors of type A/Ao containing more parameters. 
For the fermionic potential V/o hi 

V f (\,T) = x f V f0 (\)e- a ^ T2 (2.29) 

we consider two different choices. The first one is obtained directly using (2.18)-(2.20) 

12 



Vfo 



CfjyXf 
1 



+ 



729 



C 2 
C 2 

4619 £ uv 



1 + 



27 



C 2 

c? 



U + 2xf A 



(2.30) 



C 



4619 + 1714s/ - 92x 2 ] A 2 



Here one could as well use the parameter Wo which is related to C\jy by 

x f C 2 W 



C 2 
C 2 



1 



12 



For this choice the effective potential 

V cS (X) = V g (X)-x f V f o(X) 



(2.31) 



(2.32) 



has a single maximum at finite positive A = A* for all < xj < 11/2, indicating a (possible) 
infra-red fixed point. 

The second choice is obtained introducing the confinement factor (2.28) also for the 
fermionic potential, i.e., 



Vfo 



12 



C 2 
C 2 

^0 



1 + 



27 



11- 



C 



uv 



C' 1 



11 + 2s/ |A 



(2.33) 



+ t4 ( 4619 lf " 4619 + 17Uxf " 92x ') A ' (T 



y/1 + ln(l + A) 

+ A) 2 / 3 
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Now the effective potential has a maximum only at large xj . To see this concretely, consider 
again the case (2.23). The asymptotic large-A behavior of V g — xjVfo now is A 4 / 3 Vln A times 
the function 

18476 4619(1 + £ x/) 2 / 3 - 4619 + 17Ux f - 92xj 

~^43~~ 4 243(1 + 1x^2/3 ^ • (2 - 34) 

This function is positive for small Xf, negative at large Xf (< 11/2) and has a zero at 
Xf = 3.26817. Thus there is a (possible) fixed point A* only for 3.26817 < xj < 11/2. 

Let us then discuss the IR behavior of the potentials a and k which appear in the 
tachyon DBI action. For the function k we will consider the large-A asymptotics 

k(X) ~ A~ 4 / 3 . (2.35) 

A— >oo 

This is motivated by the fact that in the action the combination K(X)/b 2 has the same 
asymptotics as 1/6 2 at large A, where b s = 6A 2 / 3 is the metric factor b in the string frame. 
To ensure that the fractional exponent limit at large A does not spoil analyticity at small 
A, we replace A 4 / 3 by (1 + #A) 4 / 3 in the expression for k(A). 

More precisely, two qualitatively different, acceptable choices for the IR asymptotics 
of a (and k) were identified in [22]. These are produced by the following two choices. The 
first choice has 

and leads to tachyon growing exponentially at large r, 

r(r) ~ r e Cr (2.37) 

r— >oo 

where C is a known constant (see Appendix B) and tq parametrises the solutions. The 
second choice is given by 



K(A)= (TtW a(A) = K(A) 2^ 
and for them the leading divergence is 



1+ ( ^o + l )A + A 2 



(2.38) 



r(r) ~ CyfF^Vl, (2.39) 



r— >oo 



where the constant C is again known and now r\ parametrises the solutions. To select this 
solution, it is required that the last term in the square brackets in (2.38) grows faster than 
A 4 / 3 . 

Finally, let us summarize our choices for acceptable potentials. We always keep V g 
fixed to the expression (2.27) and choose Vjq, a, and k as follows: 

• Potentials I: We take Vjq as in equation (2.30), so that the fixed point A* exists for 
all < xj < 11/2. For a and k we use the choice of equations (2.36), which lead to 
exponentially diverging tachyon in the IR. 
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• Potentials II: We take again Vfo from equation (2.30), but use the other choice 
(2.38) for a and k. Then the tachyon diverges as r ~ y/r in the IR. 

• Potentials I*: We use now the fermionic potential Vfo of equation (2.33), which 
contains the confinement factor. Thus the extremum exists only within the interval 
3.26817 < Xf < 11/2. For a and k we use the choice of equations (2.36), which lead 
to exponentially diverging tachyon in the IR. 

• Potentials II*: We use Vfo with the confinement factor, but use the other choice 
(2.38) for a and k. Then the fixed point exist only for large Xf, and the tachyon 
diverges as r ~ y/r in the IR. 

To fully pin down the potentials, we also need to specify the value of Wq (or jCtjv) 
which is used. We choose four reference values: 

• Wo = (and constant). This is the lower bound of Wo- Actually, exactly zero Wo is 
not acceptable because the anomalous dimensions of the quark mass and the chiral 
condensate do not sum up to zero. This case is nevertheless interesting as it is the 
limit of acceptable solutions. 

• Wo = 12/11. This is the standard choice studied in [22]. 

• Wo = 24/11. For constant Wo, this is the largest possible value, for which £tjv ~~ * °° 
as xj — > 11/2. 

• Wo (and £tjv) fixed such that the free energy automatically matches with the stan- 
dard Stefan-Boltzmann result at high temperature with the correct number of degrees 
of freedom (see Eq. (2.23) and the discussion in Sec. 3.4 below). 

An ongoing work [45] studies the meson spectra in this model. As it turns out, the 
potentials I and I* admit linear "Regge" trajectories, so that the quadratic masses are 
asymptotically linear in the excitation number, m 2 ~ n, independently of the other quan- 
tum numbers. Potentials II and II*, however, have linear trajectories only in the glueball 
sector, while the other trajectories are quadratic, m 2 ~ re 2 . As linear trajectories are 
expected in QCD, this observation favors potentials I and I*. 

2.2.3 IR fixed point and the BF bound for the tachyon 

Now that the potentials are defined, one can check that they satisfy an important require- 
ment: they permit the determination of the bulk dilaton mass and, equating this with 
the Breitenlohner-Freedman (BF) instability bound, the determination of the start of the 
conformal window. Take r(r) =0 (there is no chiral symmetry breaking in the conformal 
window) and note that at small A, V g (X) — XfVfo(X) > 0. However, Vfo(X) grows faster 
and the conformal window starts at the value A* defined by the vanishing derivative 

V^K) - x f V' f0 {K) = 0. (2.40) 
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Given A* one defines an IR AdS radius 
12 



V g (\*) - XfVfo(K), Ajv > Ar- 
The tachyon mass at A* in units of £tr becomes 



f 2 

IR 



(2.41) 



2 pi 
-m m L, m 



24a(A*) 

«(A,)[y g (A„)-x / y f o(A,)] 



Gravity solutions with r = are stable when mf R Cf R 
starts when (2.42), function of Xf, has the value 4. 



(2.42) 

4; the conformal window thus 





PotI 


PotI* 


Potll 


Potll* 


W = 


4.10209 


4.33334 


4.17825 


4.38493 


w = 12/11 


3.99591 


4.33334 


4.07968 


4.38493 


Wo = 24/11 


3.71607 


4.33334 


3.80086 


4.38493 


W SB 


3.59172 


4.33334 


3.70008 


4.38493 



Table 1: The critical values x c for the various potentials. Notice that for the types I* and II*, x c 
is independent of Wq. 



Eq. (2.42) can be evaluated for the two choices of a, k above. For the choice (2.36) 
(types I and I*) the equation becomes 

3 6[1 + ^ (11 5_ 16x/) ^]4/3 

£l v [V 9 (\*)-x f V f0 (\*)} ' 1 ' J 

For the choice (2.38) (types II and II*), the x c -equation (2.42) has the form 

36[1 + 57 (115-162^* + A* 21 



4 . (2.44) 



The values of x c can then be calculated by inserting the potential V g — XfVfo and the 
chosen value for Wo in these equations. The critical values for the potentials listed above 
are given in Table 1. 

The x f -dependence of the tachyon mass for all the potential choices suggested above is 
shown in Fig. 3. The critical value x c is the rightmost point where the curve intersects the 
horizontal dashed line where the BF bound is saturated. For potentials I* and II* (solid 
magenta curves) the fixed point only exists for x* < xj < 11/2 with x* ~ 3.27. In this 
case the tachyon mass diverges as x approaches x* from above. 

From (2.43) and (2.44) one sees, using the asymptotics of the potentials (see Eq. (4.3) 
below), that — w 2 R £ 2 R ~ 1/y ln(l/x/) for type I and — m 2 R £ 2 R ~ 1/^/ for type II as 
Xf — > 0. They thus behave completely differently in this limit, for type I the mass vanishes, 
for type II it grows without bounds. In particular, for potentials I and for low Xf the 
(absolute value of the squared) tachyon mass dives below the BF bound. This means 
that the existence of a solution with a nontrivial tachyon profile and zero quark mass is 
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Figure 3: The squared tachyon mass at the IR hxed point, see Eqs.(2.43) and (2.44). Left: 
potentials I, right: potentials II. The blue curves give the masses for constant Wq. The dashed, 
solid and dotdashcd curves have Wo = 0, 12/11, and 24/11, respectively. The dotted red curves 
have Wo fixed according to the Stefan-Boltzmann normalization of the free energy in the UV. 
The solid magenta (uppermost) curves are for potentials T and IT, for which the tachyon mass is 
independent of Wo- The black dotted horizontal line marks the BF bound. 

not guaranteed [22], which means that chiral symmetry could remain intact even at low 
temperatures. However, in most of the cases, such a solution anyhow exists all the way 
down to xj = 0, and the expected picture with chiral symmetry breaking is obtained. We 
shall discuss this issue in more detail below. 

3. V-QCD at finite temperature: equations and their solution 

The V-QCD action has two kinds of vacua at finite temperature, either with identically 
vanishing tachyon or with nontrivial tachyon profile. The tachyonless black hole solutions 
can be constructed in the same way as in the Yang-Mills case [27]. Below most of the 
discussion will in principle assume the presence of the tachyon, but the construction for 
the solutions without the tachyon can be obtained simply by setting r = everywhere. 

3.1 Equations and numerical solution 

The goal now is to find numerical solutions of the Einstein's equations for the metric 
functions b(r) = e A ( r \ f(r) and the scalars A(r), r(r), satisfying 



where marks the location of the horizon. 

Due to the singular behavior of the solutions near the UV boundary (r — > 0) , it proves 
to be convenient to use A = In b as a coordinate instead of r in the numerical solution. 
Carrying out this transformation, one finds that the combination 





q(A) = e 



A dr_ 
dA 



W 



1 



(3.2) 



- 20 - 



appears naturally. This is just a rewriting of the superpotential 



W 



b 



The equations of motion then become 

q > 4 A' 2 /' q 2 

q 3 x 2 f / 

4 A' 2 /' q 2 

12 -^ + 3 7 = 7 



_ A dA 
dr 

(v g -V f Vl + fKT' 2 /q 2 ), 

Vf ) 

^l + fKT' 2 /q 2 ' 



V a 



A l' f" 

4 - — + — 
q f 

-" i i ql j- f j- x ^ lnK , 



= o, 



(3.3) 

(3.4) 
(3.5) 
(3.6) 



r.2 



4+ 4r + — + A' ^ ; |t' 3 



2/ 2 <9A 
+ 51n_^ r/2+ g 2 ainV> 



<9A 



/k <9t 



(3.7) 



A",/'A'_A' A' 2 </A' 3g 2 A^y s glfr ^ | / ^ / V f d f x r' 2 



A + / A + 4 A A 2 g A 8 / ^ dX dX 



2q 2 



1 + ^kt' 2 

q 2 



(3.8) 



Here the prime denotes differentiation with respect to A. Near the UV boundary r = 0, 

(3.9) 



A = ln6 = ln— -»• +oo. 



The range of A thus is Ah < A < +oo, where Ah is the horizon, 

f(A h ) = 0. 



(3.10) 



Numerical integration starts by solving q , X , / , r from the four first ones in terms 
of lower derivatives; the fifth equation, the equation for A, will be used as a check and 
constraint. For brevity we introduce two square root factors: 



R 1 = Jl + l^r' 2 , 



(3.11) 



and 



R2 



'12 + 



3/' q 2 



f f 

The equations to be solved numerically then are 



V n 



Yl 

Ri 



q =q 



4 + 



f 6/ 



2Vg - VfR-y - 



Ri 



X' = -^XR 2 , 



(3.12) 

(3.13) 
(3.14) 
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f" = f 



f Q 2 ( V f 



(3.15) 



+ 



6/ 

V3( 



2V g -V f R 1 



t' + — r' 3 



Yt 

Ri 

\d\K 



4 + 



R 2 + 



q* V 2/ 



J3 



7 + 5t' 3 



Xd x V f 



R 2 + 



/« 



+ T' 



/2 



(3.16) 

^7, 



In the A equation the minus branch has to be chosen as X{A) is a monotonically decreasing 
function of A. The derivatives are with respect to A. The equations are autonomous in 
the sense that there is no explicit A dependence. Numerical integration then proceeds as 
follows: 

1. Let us fix the horizon at A = Ah = — e, where e is a sufficiently small number, 
e.g., e = 10~ 6 . the values of the functions at A = 0, which is taken as the initial value of 
numerical integration, are computed by using the expansions (B.19)-(B.22) in Appendix B. 
These numbers can now be obtained by inserting the values of A^, t^, f' h . Among these the 
horizon values of the scalars, A^, 77,, remain as parameters, f' h can be given an arbitrary 
positive value, +1, say. One then finds a solution qi(A), fi(A), Xi(A), t\(A) valid from 
A = to some large upper limit A + by using NDSolve of Mathematica. The spatial 
coordinate r{A) can then, if needed, be computed by similarly integrating the differential 
equation 



r'(A) 



l q(A) 



(3.17) 



with the initial condition r(A = 00) = 0. 

2. The so obtained first-level solution fi(A) is scaled to one in the UV (A —> 00) by 
writing /2(A) = fi(A)/ fi(A + ). Simultaneously q2(A) = qi(A) / y/ f\(A + ), which is needed 
since Eq. (3.4) demands that q 2 / f be invariant. Finally, A2 = Ai, T2 = T\. 

3. The final scaling is performed to guarantee that all solutions use the same unit of 
energy or, equivalently, have the same integration constant in the integral of the definition 
(2.8) of the beta function. This implies 



A - A = ln(6) - A, 



b 2 



'4 



X(A) + 0(A 2 



(3.18) 



where Aq is the integration constant. By inserting the UV expansions of A and A from 
Appendix A, we identify Aq = In(CxjyA). We wish to scale A to one 12 , and therefore define 



Ac 



InC 



UV 



In A 



(3.19) 



and shift solutions by Aq. In practice, one implements this by determining, for a given 
numerical solution (the 0(X 2 ) term is optional), 

'b 2 6? s 



A 







lim 

A— >oo 



A-ln£ 



1 



bi 



uv 



b X 2 (A) b 2 



HboX 2 (A)) 



bl bl 



X 2 (A) 



(3.20) 



12 After this, all quantities are expressed in units of A; omitting the factor £tjv would give a unit of 
energy depending on Xf 
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and then performing the scaling 



X 3 (A) = X 2 (A + A Q ) (3.21) 

etc. for all the functions at level 2. The set q 3 (A), /3(A), A3 (A), t 3 (A), parametrised by 
the values of A^, r^, is the final numerical solution. Note that the horizon has now been 
shifted to 

A h = -A Q - e » -A ; (3.22) 
at level 2 it was defined by /2(— e) = 0. 

3.2 Physical quantities 

The set of functions q(A), f(A), X(A), t(A) (leaving out the index 3) can now be converted 
to various physical quantities: 
The temperature is 

1 P A p-Ao 

T - ~J'^ - -IMa/ (A ^- - ^^f'i-M (3-23) 

and the value of b at the horizon is 

b h = e~ Ao . (3.24) 
The quark mass m q is defined by the UV expansion of the tachyon: 

9 

r ( r ) = £ uvm(? (-lnAr) 22 ~ 4x f r (3.25) 
so that, using the relation (3.18), 

1 (b\ 



m 



lim £tjy t(A) exp 



4 ~.4"i UV 



(3.26) 



In practice, the extrapolation to A = 00 can be carried out by measuring rh q (A), as defined 
by the right hand side of Eq. (3.26), at two large values of A and then linearly extrapolating 
to A = 0: 

_ rh q (A l )X(A 2 )-m q (A 2 )\(A l ) 

X(A 2 ) - X(A l ) ■ [6 - 11] 

Linear extrapolation is chosen, because the leading neglected terms in the expansion of 
Eq. (3.26) are (up to logarithmic corrections) linear in A. 



3.3 Fixing quark mass 

The above is for fixed A^, r^. The really demanding task is to find the field configurations 
at fixed m q . For this one needs the curves Th(Xh,m q ). The quark mass is determined by 
the UV behavior of the tachyon: r(r)/£uv — m g{~ lnr) _7 °/ f ' r. To fix m q at fixed A^ we 
have to solve the equations of motion at various and find that value of which leads 
to the desired UV behavior of r(r). 
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Figure 4: Explicit bulk configurations. Left: UV large- A region. Right: Near horizon region. For 
this configuration T = 0.3839, b h = 0.7200, m q = 0.05422. 

3.3.1 Zero quark mass 

In particular, we are interested in m q = 0. This case splits in two parts: either r(r) = 
identically (chiral symmetry holds) or r(r) nonzero (chiral symmetry broken). 

If r = 0, solutions with m q = are obtained simply by setting = above. The 
solution is then controlled by the effective potential V g (X) — XfVfo(\). For classes I and 
II, this increases monotonically from A = 0, but since Vfo grows faster, the derivative 
decreases and becomes finally zero at some A = A*(xj) (see Eq. (2.40)). The extremum 
of the potential marks the location of the IR fixed point, which is screened by the horizon 
at finite temperature. Indeed, the tachyonless black holes have < A^ < A*, and for 
Xh very close to A* we obtain configurations where the dilaton is approximately constant, 
A ~ A/[ ~ A* for a long range of the coordinate before the horizon is reached in the deep 
IR. 

For classes I* and II*, the effective potential V^(A) — XfVfo(X) does not have an ex- 
tremum for xj below x* ~ 3.27. In this case the fixed point is absent, and the tachyonless 
black hole solutions are qualitatively similar to Yang-Mills (xf — > 0) [26]. In particular A^ 
can take any value. 

For non-zero r(r), the discussion of m q = configurations has to take into account 
the existence of Efimov zeroes, oscillatory behavior when approaching r = 0, which was 
discussed above in the introduction. We discuss here the standard picture which is seen 
in most cases for Xf < x c . A rough description of more complicated cases is given in 
Appendix C. The situation is summarised in Fig. 5. For large > 0, r(r) decreases 
monotonically from towards r = and ends with positive m q . We evaluate m q using 
(3.27) with two large values of A (corresponding to a small UV cutoff e in the r-coordinate). 
When Th is decreased, ultimately an (approximate) m q = configuration (ro(r, m q = 0) in 
Fig. 5) with monotonically decreasing r(r) is obtained. 

This defines the curve Tho(Xh) in Fig. 6. One finds that these solutions are possible 
only if Xh is larger than a fixed positive value, which we call A cnc j. Decreasing further, 
r(r) first develops a zero so that m q < 0. Continuing even further we find a second location 
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Figure 5: Left: Schematic presentation of the curves r/j(Aft, m q — 0), i.e., those values of the 
horizon value of the tachyon which lead to configurations with m q = 0, at that particular A/j. 
Choosing Th on the curve Th% leads to a r(r) without a linear term and with i zeros at some r. 
The zero mass solutions with vanishing tachyon live on the line < A/j < A*. The plus and minus 
signs indicate the sign of the quark mass in each region limited by the curves Tho, Tuiy ■ .Right: 
Schematic presentation of the r dependence of the bulk tachyon for low node numbers. Tachyon 
solutions for top to bottom are: a generic solution with m q > and no nodes ("0 node"), the 
standard solution with zero quark mass (to), the solution with zero quark mass and one node (ti), 
and a generic solution with m q < and one node ("1 node"). See the text for a more detailed 
explanation. 

where m q = vanishes. This is a configuration with one tachyon node (ri(r, m q = 0) in 
Fig. 5). The pattern continues with an ever increasing number of nodes, until one ends up 
with the curve Th c (Xh), below which a solution with the standard UV boundary does not 
exist. Numerically, the curves r/ l o(A/ l ) and Thi(Xh) can be separated, but already r/ l 2(A/ l ) 
would require so much effort that we have not embarked on computing it. As we approach 
the conformal window, the curves t^q, t^i, ■ ■ ■ get closer and closer to 77^ and finally vanish 
for Xf > x c . 

We expect that increasing the number of nodes increases the free energy so that to 
study equilibrium states it is enough to compute r/ l o(Aft,m (? = 0). This was checked at 
zero temperature in [22] numerically for potentials I, and analytically in the limit Xf — > x c 
as well as in the limit of large number of tachyon nodes. 

3.3.2 Nonzero quark mass 

For nonzero quark mass the special solution with identically vanishing tachyon profile is 
missing. However, there are solutions of various types for r/j > 0, as suggested by Fig. 5. 
We shall here restrict to the "standard" solutions which have monotonic tachyon, i.e., the 
region above Tho in Fig. 5 (left). Below this curve there can be Efimov type solutions 
where the tachyon has nodes. As for m q = 0, we expect that these solutions have higher 
free energies than the standard one. In the region of standard solutions, the dependence 
of quark mass is smooth (see Fig. 6). We have found numerically that for fixed Xh the 
correspondence between m q and is one-to-one. Therefore m q can be kept fixed by 
following a set of well-defined curves on the (A^, T/J-plane, some of which are sketched in 
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Figure 6: The curves Th(Xh,m q ) for values of m q marked in the figure, computed for Potential I 
with W = 12/11. Here A end = 0.5221 and A* = 0.6467. 

Fig. 6. 

It is also interesting to notice how the m q = solution is obtained from the ones having 
finite quark masses as m q — > 0. What happens for nonzero quark mass is shown in Fig. 6 for 
a concrete computation. If < < A cn d (and A^ fixed), the curve Th(Xh,m q ) approaches 
zero as m q — > 0, indicating that r(r) approaches the chiral symmetry conserving solution 
(r(r) = 0) uniformly. If A^ > A cn d, ^(A^,™,) approaches Tho(Xh) instead, which implies 
that r(r) converges to the standard chiral symmetry breaking solution ro(r). 

3.4 Thermodynamics 

We now want to compute minus free energy density or pressure p(T, m q ; xj) of the gravity 
dual, assuming that all the quarks have the same mass m q . In particular, we are inter- 
ested in m q = 0. The chemical potential is zero, there is an equal number of quarks and 
antiquarks. The equilibrium phase has the largest pressure. 

The basic strategy is to compute the temperature and entropy density from the for- 
mulas 

T= "ir //(r/l) ' S= lGf {rh) ' (3 - 28) 

where / and b are obtained by solving Einstein's equations. The pressure is then obtained 
by integrating s(T) = p'(T). The key technical issues are keeping track of the quark mass 
and specifying the integration constant in the pressure integral. 

The general structure of temperature (for a case containing a fixed point) is shown in 
Fig. 7, to be consulted in association with Figs. 5 and 6. For m q = two branches separate. 
Firstly, for < A^ < A* there is the temperature computed for chirally symmetric vanishing 
tachyon solutions. We shall use the notation T u (Xh) = T(A/ l ,r/ l = 0) for this temperature 
below. 

The chiral symmetry breaking solution exists for A cn d < Xh < oo and as X^ — > A en d, 
the corresponding temperature curve ends precisely on the curve which has identically 
vanishing tachyon. The temperature curve is computed by using the zero node zero mass 
curve Tho{Xh) in Fig. 5. We shall use the notation ^(A^) = T(A^, t^A/j, m q = 0)) for this 
temperature. If we computed the temperature for the one node solution r^ 1 (A/ l ), we would 
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Figure 7: The temperature as a function of A^ for solutions for Pot II at Xf = 3 and Wo — 12/11, 
both for zero (left) and very small mass (right). The asymptotic limits (3.29) and (3.30) are also 
shown for m q = 10 -5 , in the range of the figure the UV limit is not yet accurate. The maximum 
value A* of A for the r = curve is defined in (2.40). See also Fig. 2. 



get a curve which lies significantly below the zero node curve in Fig. 7 and again ends on 
the zero tachyon curve. These solutions will have a higher free energy and we can thus 
neglect them. 

Whenever the quark mass is nonzero, the tachyon cannot be vanishing and that branch 
disappears. However, as seen from Figs. 6 and 7, the small-mass curve very closely approx- 
imates the zero tachyon curve, also at small A^. 

Analytic approximations are often useful. In the UV f(r) — > 1 — ?" 4 / r h so that 

*T(Afc) = - = e l ^(b X h ) b ^ = (3.29) 

Similarly, in the IR (see (B.10) in Appendix B), 

r(A fc )~(|lnA ft ) 1 /2 j b (\ h )~ -1^(2 In Aft) 1 /*. (3.30) 

\ 

For a numerical check, see Fig. 7. The interesting physics takes place in the region con- 
necting these two limits. 

The function b(Xh) decreases monotonically while the function T(Xh) decreases in the 
UV but starts increasing in the IR. The physics of the UV increase is obvious, this is 
the weak coupling limit which naturally corresponds to large T of a thermal fluid. The 
(extremely slow) increase in the IR is a quantitative fact but does not correspond to a 
stable phase. This is simplest seen by computing the sound velocity 

2 dp s b(Xh) ( dT 



° s de Ts'(T) 3T(-V(X h ))\ dX J ' (3 ' 31) 

A stable phase has c 2 s > (equivalently, has a positive specific heat) and this requires 
T'(Xh) < 0. Thus only the UV decreasing part can correspond to a stable phase, the 
IR part is the unstable small black hole region, small since s ~ b 3 (Xf t ) — > there. It is, 
nevertheless, crucially important for the phase structure. 
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To compute the pressure, we have to integrate the entropy density (3.28) over T. 
Taking Xh as a variable, we have integrals over the two branches in Fig. 7: 

1 f°° 

Pb( T ) = J7T I dX h (-Tl(X h ))bl(\ h ) +p b (oo), (3.32) 

3, 



4^5 Jx h (T) 



1 



A,. 



Pu(T) = / d\ h {-T' u {\ h ))bi{\ h ) + p u (A*), (3.33) 

4< -5 Jx h (T) 

where b, u refer to the chiral symmetry broken (r^ = Tho(Xh,m q = 0)) and chirally sym- 
metric (or unbroken, r = 0) phases. The continuity of pressure at T en( \ = T(A en d) leads to 
a rather remarkable consistency check of the entire scheme: it demands 

— / dx h {-r u {x h ))bl{x h )-— / dx h (-r b (x h ))bl(x h )= Pb (oo)- Pu (x*). (3.34) 

However, the difference on the RHS is nothing but the difference between the free energies 
of the broken and symmetric phases at T = 0: 

Pb (oo) - Pu (X«) = -F b (T = 0) + F U (T = 0). (3.35) 

This difference was computed in [22] from the T = solutions, with no black hole. Here 
they are computed in (3.34) from the black hole solutions and we have checked numerically 
that the results agree within the numerical precision. 

The computation of the free energy now proceeds as follows, first for the simple struc- 
ture of T(Xh) in Fig. 7: 

• Start by integrating (3.32) from some large value of down to A cn d, choosing 
Pb(co) = 0. Since T'(Xh) > in Fig. 7, this leads to a negative pressure. This 
is not the stable phase, the physical stable phase is not described by this metric. The 
stable phase with the largest pressure is the thermal gas phase with p = 0. 

• At A cn d move to the chirally symmetric r = branch and fix the constant ^(A*) by 
demanding continuity of pressure. Since now T'(Xh) < 0, p starts increasing. At first 
p is still negative and the stable phase is the thermal gas phase with p = 0. 

• At some A^ = A c pressure passes through 0. This defines a transition temperature 
T/j since from now on the black hole metric has the largest pressure. Since r = this 
black hole phase is chirally symmetric. 

• The latent heat of the transition is 

L s(T h ) 1 /6(A^)\ 3 2 4^ 2 7 
^ = ^T = 4G5V^W <N ^ (l + * Xf) > (3 - 36) 
where the maximum value is obtained taking normalisation from (3.39) and using 
the UV approximation (A. 19). Counting degrees of freedom one has Nj Goldstone 
bosons in the low T phase (for which we do not have a T dependent gravity dual) 
and 2N% + ^ N c Nf degrees of freedom in the high T phase. These are equal at xj = 4 
and if latent heat is naively assumed to be proportional to the jump in the number 
of degrees of freedom, one might rather expect L to decrease when xj increases. 
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Asymptotically, for large T, — > we have ttT = \jr\ l = b(Xh) / £tjv so that 

/•oo 

4G 5P = (nCwf / dx(-T'(x)) T 3 (x) = \ (vr£ uv ) 3 r 4 . (3.37) 

If one for large T assumes that the system becomes a gas of non-interacting bosons 
and fermions one should have 



This is obtained from (3.33) if 

1 A. 1 -I- I T. 4 

-iV c 2 , (3.39) 



4G 5 45vr £3 



'UV 

which can be used to normalise the pressure. 

• The above was for the simple T(A/J in Fig. 7. Depending on the potentials, more 
complex structures can appear, as analysed in the following section. 

• To present results for p/T we choose to normalise it so that it approaches at large 
T the ideal gas Stefan-Boltzmann pressure according to (3.38). However, we have no 
dynamical argument for fixing the xj dependence of Cjjv m (3.39). We shall present 
the phase diagrams for two choices, for the automatically SB-normalised case (see 
Eq. (3.39)) 



4 ' 



£uv = (l + Ix / ) 1 / 3 , W = — 

Xf 



1 



(1 + fx,) 2 / 3 



(3.40) 



and for the Wq fixed case 



1 12 24 

Ajv = , W o = 0, — ,— (3.41) 



In the former case one simply has 

1 



4G 5 45vr 

and in the latter case 13 



Nl (3.42) 



1 + lXf AT 2 ; (3.43) 



4G 5 45vr (1 - i Xf W ) 2 / 3 
the factor iV 2 is furthermore often implied, i.e., results for p/(N 2 T ) are given 



13 Notice that in this case the glue part of the V-QCD action will also depend oni/ through the normal- 
ization factor 1/4G?5. 



- 29 - 



SU(N) gauge theory, massless fermions 

T 















Conformal 






\ 


Window 




X broken 


\ 

\ 

\ 


(X symm.) 






\ 








\ 
\ 








X 

\ 








\ 














1 


2 3 J 


i c 4 5 X 


IS 



-V 



Figure 8: Qualitative behavior of the transition temperature between the low and high T phases 
of V-QCD matter. 

4. Results for the phase structure 
4.1 Phase transitions 

Let us first review what one qualitatively expects for the phase structure of V-QCD when 
the number of (massless) fermions is changed [46]. This is shown in Fig. 8, where the 
transition temperature between a low T and a high T phase is plotted as a function of Xf. 

A few reminders are in order. In the absence of quarks, YM has a Z^ c center symmetry 
that is central in the definition of the confined and deconfined phases. The relevant order 
parameter is the Polyakov loop that transforms nontrivially under Z^ c . If its expectation 
value is zero, we are in the confined phase, while the expectation value becomes non-zero 
in the deconfined phase. 

This expectation value is simple to calculate holographically, [47]. It corresponds to 
a string world-sheet along the time circle, and hanging down straight in the holographic 
(radial) direction. The important difference is where it ends. At zero temperature, this 
worldsheet extends to r — > oo and is the world-sheet of a free quark. Standard renormal- 
ization subtracts its contribution completely and therefore the Polyakov loop vev is zero 
(to leading order in 1/N C ) in the zero temperature phase. 

In a regular black-hole phase, the worldsheet terminates at the horizon and after 
subtraction the Polyakov loop expectation value is non-zero. This is in agreement with the 
identification of black- hole phases generically as deconfined phases. 

In the presence of massless quarks, the center symmetry is not a symmetry any more, 
and the Polyakov loop is not an order parameter. However at large iV c , there is alternative 
order parameter for a deconfined phase, namely the N c dependence of the free energy, F. 
In the confined phases F ~ 0(1) while in deconfined phases, F ~ O(N^). Again, with 
this criterion, the vacuum solutions (without horizons) are "confining" (F ~ 0(1)) while 
any black hole solution with regular horizon is "deconfined" (F ~ O(N^)). It is therefore 
natural to use this criterion in our analysis in order to define deconfined phases. 

The true symmetry in the case of massless quarks is chiral symmetry. This always 
has an order parameter, the chiral condensate, that distinguishes chirally symmetric from 
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Figure 9: Examples of the Tend: 2fc an d 2c roS sover transitions in potential II with Stefan-Boltzmann 
-normalization of £\jy and with xt = 3 (see also Fig. 13). Upper: The temperature T(A^) . The 
curving of T u (Xh) at A^ ~ 0.2, T ~ 2 is related to the crossover. The inset shows the minimum 
of Tb(A/j), which causes p& to be positive between Th and r en( j. For comparison, we also plot 
T(Xh) for IHQCD with Xf = 0. Lower left: p/T 4 in a close-up around the region of the T/j and 
^end -transitions. Lower right: an overview of the pressure in the same case, also showing the 
interaction measure, the peak of which determines the position of T crossover . The black curve shows 
the vacuum beta function, scaled to fit, as a function of temperature in the symmetric phase, so 
that j3(T) = /3(A U (T)), where X U (T) is the inverse function of T u (Xh)- The walking maximum 
of the beta function clearly coincides with the plateau related to T crossover , confirming that the 
p/T ~ constant phase below T crossover is indeed the quasi-conformal phase related to walking 
dynamics. 

chirally broken phases. 

Given the remarks above, we summarize what we would expect. 

• For Xf = one has the Yang-Mills 1st order phase transition between a confined and 
deconfined phase. In the high T deconfined phase, the Z^ c symmetry is broken. 

• For a somewhat higher Xf one expects that there still is a 1st order transition. How- 
ever, now this transition will involve chiral symmetry breaking/restoration. 

• For x j approaching x c one expects the transition temperature to decrease rapidly as 
follows from Miransky scaling. 

• For Xf in the conformal window, x c < Xf < 11/2, both the low and high T phases are 
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Figure 10: An example of the T s transition in potential I with Wq = 24/11 and with Xf = 0.3 Left: 
The local maximum and minimum which generate the 1st order T s -transition. Right: p(T)/T in 
the region around which the 1st order T s transition takes place, extending to smaller T in order to 
show the relation to the and T en( ^ transitions. 

conformal ones, which can be separated by a crossover. The only transition happens 
at T = + like in the AdS black hole in Poincare coordinates. 

The models we consider contain the full fermion backreaction and therefore predict a 
somewhat more detailed phase structure. New phase transitions of different orders can take 
place, lines can split in two, etc. The behavior in the conformal window (x c < Xf < 11/2) 
is nonetheless always simple: there are no transitions, but a crossover between the low and 
high temperature conformal phases. Therefore we concentrate first on the phase structure 
in the region below the conformal transition (xf < x c ). 

While the details of the phase structure depend on the choice of potential, the various 
phase transitions encountered appear in certain systematic ways. We will define a consistent 
notation, and describe the classes of transitions, assuming the system is heated up and we 
go from low temperatures to high temperatures. 

To motivate the notation, we first list the various transitions and the corresponding 
temperatures. 

• Th is the analogue of the QCD hadronisation transition if it is the chiral restoration 
transition (chirally symmetric — > chirally broken). 

• T en d is the end point of the curve T^A^) = T(Afe, t/jo(A^, m g = 0)), which contains 
the black holes with tachyon hair. For values of A^ smaller than at this endpoint, the 
black-holes have no tachyon hair. 

• ^crossover marks the position of a crossover. This crossover is defined by the position 
of the peak in the equation-of-state ((e — 3p)/T 4 ) as a function of temperature. 

• T s takes place at small Xf within the chirally symmetric phase when one can jump 
from one decreasing branch of T u (\h) (no tachyon hair) to another. 

• Finally Tyi involves the splitting of one 1st order line to two. 



- 32 - 



With this notation we may now describe in detail the various types of transitions and 
crossovers we have found, and show examples of each case. In the figures we denote the 
stable phases with thick lines and meta- and unstable phases with thin lines. 

• The 1st order hadronisation transition at T^, happens either between the chirally 
broken — > a chirally symmetric phase (see Fig. 7) or from a chirally broken — > a 
chirally broken phases (see Fig. 9). 14 As described above, our normalization for 
pressure is such that the pressure of the (T = 0) hadron gas phase is zero. In the 
holographic setup, this transition is between that of the black hole phases, whose 
pressure remains positive down to the lowest temperature, and the hadron gas phase. 
The transition takes place at the temperature Th where the pressure of the BH phase 
reaches zero. Whether this phase is chirally symmetric or non-symmetric depends on 
the potential choices and Xf. For an example, see Fig. 13. 

• The 2nd order chirally broken — > chirally symmetric transition at T en( j = T(A en d ; 0) = 
T(\ en d, Tft,o(A en d) "iq = 0)), see Figs. 7 or 9. Since the chiral symmetry breaking 
solution starts to exist only above some A en d> the system makes at that point a 
transition to the chirally symmetric phase. However, this transition may be absent in 
the thermodynamic limit: if Pb(^h) is everywhere negative, the transition is between 
two thermodynamically metastable phases, and the relevant saddle point is never 
dominant. We denote the temperature of the transition by T enc q. Since this transition 
takes place at one single value A/j = A enc j, both pressure and entropy density are 
continuous (b(Xh) does not jump). Therefore, only p"(T) or (? s are discontinuous, 
and the transition is of second order. 

• The high-T chirally symmetric — > chirally symmetric crossover at T crossover , see Fig. 9. 
This is a crossover which is expected on general grounds when xj is near but below 
x c . It reflects the change of the dynamics from the walking region, where the QCD 
coupling constant evolves slowly, to the region in the deep UV where it runs. In this 
sense, above the crossover it is the nontrivial fixed point theory that controls the 
thermodynamics, while below the crossover it is the YM-like theory that controls the 
dynamics. 

The thermodynamics behaves as follows: At first p/T 4 stabilizes to some intermediate 
value, before eventually increasing very slowly toward the Stefan-Boltzmann limit. 
For the potentials studied here, this creates a clear, although very broad, peak in 
the interaction measure, and the position of that peak can be used to define the 
temperature T crossovcr at which there is a crossover. The peak of the interaction 
measure is also observed at low values of Xf. In this region, however, T crossovcr is 
typically relatively close to Th- Note also that for SU(iV c ) YM theory, Nf = the 
interaction measure starts decreasing immediately at Th [48], T crossovcr ~ T; L . 

14 There is also the special case of potentials I» at low Xf where the transition analogous to Th takes place 
from a chirally symmetric thermal gas to chirally symmetric black hole phase (see Fig. 19). 
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Figure 11: An example of the T\i transition in potential I with Wq = 12/11 and with Xf — 3.5 
Left: The overall structure of T(Xh), with an inset showing the maximum and minimum in more 
detail. Right: a close-up of p(T)/T 4 in the region where the T12 -transition happens, with an inset 
showing further detail. 

• The 1st order high-T chirally symmetric — > chirally symmetric transition at T s , see 
Fig. 10: With some choices of potential, at low xj, T U (A/J in the chirally symmetric 
(unbroken) part of the solution develops a local maximum and minimum. There are 
then two values of Xh between which both the pressure and the temperature of the 
solution match, and there is a 1st order transition between these two branches of the 
chirally symmetric solution. Interestingly, T s approaches the temperature of the YM 
transition in IHQCD as Xf — > (see the discussion in Section 4.8). 

• The 1st order chirally broken — > chirally broken transition at T12, see Fig. 11. This 
happens in the chirally non-symmetric phase, with potential I and Wo = 12/11, T(Xh) 
which develops a local minimum and maximum at large Xf. This again induces a 1st 
order transition, which we denote by T12. In this case the single 1st order transition 
at Th splits into two 1st order transitions as ij increases above some critical value. 
Above this value, the transition with higher (lower) temperature is identified as T12 
{T h ). 

4.2 Class-II Potentials 

Let us then discuss the details of the phase structure for the various potentials and choices 
of W defined in Sec. 2.2.2. 

We take Class-II first since it leads systematically to a simple phase structure. We 
observe two possibilities: First, for xj up to some value x x < x c the 1st order deconfinement 
and chiral transition temperatures coincide, Tj, = T x , from this value up to x c one has 
T x > Td and the higher chiral transition is of 2nd order. Second, Td = T x all the way up 
to x c and x x is absent. 

For this choice of potentials the tachyon diverges ~ y/r — r\ at large r. The part 
Vfo(X) of the fermionic potential is given by Eq. (2.30) and a(A) and n(X) are given in 
(2.38). Notice that the deconfinement temperature always equals the temperature of 
the "standard" 1st order transition in the holographic framework. The temperature of 
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Figure 12: An example of a configuration where all but the crossover and hadronisation transitions 
^crossover, Th, are in the thermodynamically unstable region, in the initial stages of the approach 
to the IHQCD limit. The potential is II with Wo — 12/11 and with Xf — 0.4 Le/£:The temperature 
T(A/i). Note that everything to the right of the Th transition is in the unstable phase. Right: 
p(T)/T A in the region where the Th transition and the unstable T en( ^ and T s -transitions happen. 

the chiral symmetry restoration T x can be either T en d or Th depending on the order of the 
transitions, see examples below. 

The result for the SB-normalised case is shown in Fig. 13. For < Xf < x x ~ 2.46 
we find that T en( j < Th, but r enc j is in the metastable branch of the solution. Thus the 
deconfinement and chirality transitions coincide here, T^ = T x = Th- In other words, if one 
could sufficiently supercool the system below Th in the high-T chirally symmetric phase, 
the symmetry breaking transition could take place at T em j < Th- In the thermodynamic 
limit there is no supercooling and only Th is seen. 

Above Xf ~ 2.46, the second order T en d moves above Th and becomes stable, as seen 
in the bottom right plot of Fig. 13. Therefore, we first have a 1st order Th transition from 
the thermal gas solution to a chirally breaking black-hole phase, and then a 2nd order 
transition from the chirally broken low-T phase to the chirally symmetric high-T phase. In 
other words, T x > T^ with a 2nd order chiral and 1st order deconfinement transition. For a 
more detailed view of the thermodynamics in this region at Xf = 3, the reader is guided to 
the left panel of Fig. 24 where the chiral condensate as well as the energy and the pressure 
are plotted as functions of T. The chirally symmetric crossover transition T crossover is for 
all Xf, the highest temperature transition. 

For x f — > x c both T en d and Th are expected to approach zero as specified by Miransky 
scaling. Numerical results are compatible with this. 

When Xf — > one would expect that the Th transition smoothly approaches the 
transition temperature of large N c hot Yang-Mills theory. Note, however, that strictly 
speaking the limit of YM theory demands Nf = and falls outside the Veneziano limit 
Nf — > oo of QCD. Thus it is not surprising that nontrivial metastable structures appear 
at Xf — > 0. What happens is that the curve T = T u {\h) of the chirally symmetric phase 
suddenly at Xf ~ 0.2 develops a local minimum similar to the one shown in red in Fig. 10. 
Further evolution of this minimum is shown in Fig. 22. Associated with this there is a first 
order T s transition in the metastable branch. It is so slightly below Th that it is not visibly 
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separated in the bottom left plot of Fig. 13. As discussed in section 4.8, both Th and T s 
approach the transition temperature of YM as Xf — > 0. T en( j crosses above all of the other 
transitions for low Xf, but it is also in the metastable branch, see Fig. 12 for details. 

The phase diagram for potential II at Wq = 24/11 is shown in Fig. 14. The phase 
structure is qualitatively similar to the SB-normalized case. For Xf < x x ~ 3.19 the stable 
Th transition is the only one in the thermodynamic limit, with T en( ^ < Th in the metastable 
branch of the solution. Thus again = T x . Above xt — 3.19, the second order T en d moves 
above Th and becomes stable, see bottom right plot of Fig. 14. Thus we again have T x > T c i 
with a 2nd order chiral and 1st order deconfinement transition. The chirally symmetric 
crossover transition T crossover is for all Xf the highest temperature stable transition, except 
between x/~ 1 to xt ~ 2.7, where the interaction measure does not have a maximum and 
the crossover therefore does not exist. 

Now T s which appears in the metastable branch slightly below Th in Fig. 10 (bottom- 
left) visibly separates from Th- Again T s and Th approach the temperature of the YM- 
transition in the Xf — > -limit, as discussed in section 4.8. T en( j crosses above the Th 
transition for Xf < 0.34, but it is also in the metastable branch, see Fig. 12 for details. 

The phase diagram for potential II at Wo = 12/11 is shown in Fig. 15. The main 
difference with respect to the previous cases is that T cnc j < Th for all values of Xf, so 
the region with T x > Tj does not exist. Notice that T en d is close to Th for Xf — > x c as 
seen from Fig. 15 (left). Because the region with small x c — xj is numerically challenging, 
we do not have reliable data for Xf > 3.8. However, nontrivial structure apart from the 
Miransky scaling, such as rapid changes in the ratios of the various temperatures, are not 
expected in this region (see discussion below in Sec. 4.8). The chirally symmetric crossover 
transition T crossover is the highest temperature stable transition where it exists. The next 
stable transition is everywhere Th, and as already pointed out, T en( j is in the metastable 
branch of the solution. Details of further metastable structure at small x f are shown in 
the right hand plot. At Xf ~ 0.25, the first order T s transition appears in the metastable 
branch slightly below Th, see Fig. 10. This transition develops into the YM transition in 
the Xf -limit. T en d crosses above the Th transition, but it is also in the metastable 
branch, see Fig. 12 for details. 

The phase diagram for potential II at Wq = is shown in Fig. 16. For all points shown, 
T en d is below Th and in the metastable branch. The crossover exists when Xf > 3.6 and 
again between Xf = to ~ 0.7. The close-up of the small Xf -region in the right hand plot 
shows the crossover and the hadronisation transition Th, with the T en( j and T s transitions 
in the metastable branch. As a new feature the crossover also becomes metastable for 
0-5 < x f < 0.7. 

Finally, let us comment on the Xf dependence of the transition temperature(s). For 
SB normalised Wo or Wo = 24/11 (Figs. 13 and 14), Th and T en d decrease with increasing 
Xf, in qualitative agreement with estimates based on field theory [49]. Decreasing Wo to 
12/11 (Fig. 15), however, the Xf dependence becomes almost flat, and for Wo = (Fig. 16) 
the temperatures increase with Xf up to Xf ~ 3.5. Rather similar behavior with varying 
Wo will be found for potentials I below, where the x j-dependence is partially disturbed by 
the additional structure appearing at low Xf. 
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Figure 13: Upper: The phase diagram for potential II, Wo Stcfan-Boltzmann normalized, x c — 
3.70. The dashed boxes show the regions detailed in the bottom two plots. In the bottom left plot 
T s Th at Xf ^ 0.2 is not visibly separated. For discussion, see text. 




Figure 14: Upper: The phase diagram for potential II, Wo — 24/11, x c — 3.80. The dashed boxes 
show the regions detailed in the bottom two plots. For discussion, see text. 




Figure 15: Left: The phase diagram for potential II at Wq = 12/11, x c — 4.08. The dashed box 
shows the region detailed in the other plot. For discussion, see text. 
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Figure 16: Left: The phase diagram for potential II at Wq = 0, x c = 4.18. The dashed box 
shows the region detailed in the other plot. For discussion, see text. Tcrossover continues into the 
conformal window. 



4.3 Class-II* Potentials 

In this section, we consider the phase diagram corresponding to the potential II*. Recall 
that the star subscript refers to the fact that the potential V(A) = V g {X) — xjVo/(A) has 
an extremum only for x* < Xf, while for the cases discussed earlier such extremum exists 
for all values of xf, see Sec. 2.2.2 for detailed definitions. 

The resulting (xf,T) -phase diagram is shown in Fig. 17, the top panel shows how the 
phase diagram is derived at Xf = 2.5. Starting at large T one is in the tachyonless black 
hole phase (thick red curve). At Th ~ 0.8A pressure goes to zero and the ground state is 
the thermal gas phase with p = 0. If one could supercool further one would at T enc j meet 
the chirally broken tachyonic black hole phase. It has a higher free energy than the stable 
broken phase and therefore is unstable. 

The main features are that the crossover exists only for small values of Xf, Xf < 2 
where it nearly coincides with Th, and again at larger values xj > 3.5, where it is clearly 
separated from Th- The second order endpoint T en( j remains in the unstable phase for 
Xf < x c . Below the conformal window, for values 2 < xj < 4 both Th and T en d increase. 
They reach their maximum and finally start to decrease (as predicted by the Miransky 
scaling) only around xj = 4, very near the boundary of the conformal window. This 
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Figure 17: Phase diagram for potential II* with Wo SB normalised (bottom). The top panel 
shows T = T(Xh) and p/T 4 at Xf = 2.5. T crossovcr reappears at about xt = 4 at a temperature 
well above the range shown here, and continues into the conformal window. 



suggests that the modification of the potential has the tendency to "squeeze" the walking 
region. 

4.4 Class-I Potentials 

For class I potentials Fig. 18 shows phase diagrams for Wq = 0, 12/11, 24/11 and for the 
SB-normalised case. Recall that for these potentials the tachyon diverges exponentially in 
the IR. The choices of a and k are given in Eqs. (2.36). We also remind that transitions 
between stable phases are plotted as thick lines. Transitions plotted as thin lines can be 
seen only if the system is, e.g., supercooled, so that they are not there in the thermodynamic 
limit. 

One can observe several characteristic features for varying Wq: 

• The first observation is the striking structure near Xf = which is observed at large 
Wo, i.e., for Wo = 24/11 or SB normalized. The temperatures and T en( } drop 
rapidly with decreasing Xf near Xf = and reach zero at a finite value of Xf. Below 
this critical value, all phases are chirally symmetric. 

This behavior is related to the tachyon mass at the IR fixed point, shown in Fig. 3. 
For PotI (the absolute value of) the squared tachyon mass is below the BF bound for 
low values of x / . Therefore it is not guaranteed that a solution with zero quark mass 
and nontrivial tachyon profile exists (at any temperature) in this region. For large Wo 
it actually turns out that the solution with m q = and nontrivial tachyon profile does 
not exist for very low Xf, which explains the absence of chiral symmetry breaking. 
This implies that this potential is not describing a QCD-like theory. However, the 
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applicability of PotI can be rescued by a simple logarithmic modification of n(X), see 
Section 4.6 and Fig. 20. 

• The symmetric—)- symmetric transition T s becomes a stable transition when Wq = 
24/11 or SB normalized. For comparison, for Potll it was always in a metastable 
phase. This happens mostly in the region of very low Xf where all phases are chirally 
symmetric so that and T en( j are absent. For Wq = 24/11 we observe a region with 
0.25 < Xf < 0.45 where these transitions are also present. In this case the order 
of transitions is T s > T en d > Th, and chiral symmetry is broken in the middle one. 
For Wq SB normalized we find instead a region with 0.2 < xj < 0.5 where only the 
crossover exists, so that the phase structure is similar to the conformal window. 

• At large Xf ^ 3, Wq = 12/11, one observes the splitting of the 1st order line Th into 
two 1st order lines T\i > Th,. The order of the transitions is T cnc i > T\2 > Th, chiral 
symmetry is broken at the largest one, T cn( j. The holographic action therefore gives 
two consecutive 1st order transitions within the chirally-broken phase. It is an open 
issue what the nature of these transitions is. It is plausible that PotI at large Wq is 
not related to QCD-like theories. 

• The high temperature crossover exists over a larger and larger range when Wo in- 
creases and ultimately appears at all Xf. This is the same tendency seen also for 
potentials in the II class. 

4.5 Class-I* Potentials 

Finally, we present the phase diagram corresponding to the potential I* in Fig. 19. The 
striking difference between the phase diagram of the potential I* in comparison with po- 
tential II* considered earlier is that for small values of Xf < 2 there are no solutions with 
broken chiral symmetry, not even at low temperatures; all phase boundaries here are be- 
tween chirally symmetric phases. There is Th, but now it describes a chirally symmetric 
— > symmetric transition. To illustrate this we show explicitly T = T(Xh) at xt = 1. It is 
very structureless, and has no solutions with nonzero tachyon. Thus the (A/,T) -diagram 
is qualitatively similar to the Yang-Mills case [27]. Only above xj ~ 2 and below x c is 
chiral symmetry broken at low temperatures. 

Otherwise the overall features are similar to those in the case of potential II*. For 
small values of Xf, xj < 2, the crossover nearly coincides with Th . The second order 
endpoint, T en( j, is in the unstable branch for small values of Xf, but enters into the stable 
branch at x ~ 4. Below the conformal window, for values 2 < Xf < 4 both Th and T en d 
increase. They start to decrease only at Xf ~ 4, very near the boundary of the conformal 
window. 

We have also studied the potentials I* for the case of fixed Wq and found qualitatively 
similar results for the phase structure for Wq = 12/11, 24/11. For Wq = the problematic 
region without chiral symmetry breaking is absent, and the phase diagram is similar to 
Potll*. This implies that, like Pot I, this type of potential is probably not applicable for 
QCD-like theories when Wq is large. 
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Figure 18: Phase diagrams for potential I. Top: Wo = 0, Middle ones: W = 12/11, 24/11, 
Bottom: SB normalisation of Wq. A blow-up of the small Xf region is shown at right separately 
for three top rows. The leftmost vertical line in the bottom three figures denotes the value of Xf 
below which chiral symmetry breaking solutions do not exist. 



4.6 PotI with logarithmic correction to k(A) 



The function k(A) in the action (2.2) represents the effects of going from the string frame 
(to which the derivation of the DBI action as the a' — > limit of open strings leads) 
and the Einstein frame (where the gravity dual is formulated). Extending the conformal 
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transformation relating these to UV by A — > 1+A one has, in terms of the metric functions, 15 

where b and b s are the metric factors in the Einstein and string frames, respectively. 

The potential (2.27) carries the factor A 4 / 3 , but also the logarithmic factor (In A) 1 / 2 , 
which plays a quantitatively important role: for (In X) p the excitation spectrum is m ~ n p 
[26] and one wants the Regge-like spectrum, P = 1/2. Also numerically lnA-effects are 
important, see Fig. 7. To study these effects in k we use the parametrisation 

[l+lnp+A)]" . . 

" (A) = [1 + 1(160 + 1+^] V (4 ' 2) 

There are constraints on this parametrisation from the UV and IR. First, to maintain 
the proper mass anomalous dimension equation (2.26) at small A, p has to appear also in 
the denominator as shown in (4.2). Secondly, for p = the tachyon grows exponentially 
in r according to Eq. (2.37). The effect of p on this comes from the change b 2 /n = 
Vln A/(ln X)P- ~ r 1-2 ' 1 (in the IR r ~ y/\n A, see (B.4)). This effect propagates through the 
computation of the r dependence which comes out to be r(r) ~ exp(Cr 1-2 ^), indicating 
that p < 1/2. 

The most interesting effect comes from evaluating the tachyon IR mass using (2.42). 
The result is shown in Fig. 20, to be compared with Fig. 3. The difficulty with PotI was that 
at small Xf the curve in the left panel of Fig. 3 dropped below the BF bound. The reason 
for this is easy to see analytically by studying the A* — > oo limit of (2.42), which gives 
— m lR^lR ~ ( — hixj) -1 / 2 in this case. For small Xf, A* approaches infinity and obviously 
negative values of p increase the tachyon mass — w, 2 R ^ 2 R , so that — "i 2 R ^ 2 R ~ (— In xj)~^ _1 / 2 
For p < — 1/2 it even grows without bounds as for Potll in Fig. 3. This is seen in Fig. 20. 

As a consequence, the phase diagram for PotI with log-modified k does not suffer 
from the problems at small Xf described earlier for PotI. The phase diagram computed 



15 Notice that we introduced additional constants in the formulas (2.36) and (2.38) in order to match with 
the perturbative anomalous dimensions in QCD. 
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Figure 20: Left: The tachyon mass at A* for PotI with log-modified k(X), evaluated using Eqs. 
(2.42) and (4.2) for fx = 0, — 4 , — h . The fi = curve is the same as the red dotted curve in left 
panel of Fig. 3. Right: The phase diagram for PotI with SB-normalised Wo with k(X) given by 
(4.2) with fi = — |. T en( j intersects at Xf = x x = 0.72. T crossover is the same as in Fig. 18, 
bottom. 

for p, = — ^ is shown in Fig. 20 and, in fact, resembles qualitatively those for Potll. This 
is very gratifying since PotI also leads to a Regge-like particle spectrum [45]. PotI with 
log-modified k(A) (4.2) thus seems to be the gravity dual leading to the simplest thermo- 
dynamics in Fig. 8 and the expected Regge-like hadron spectrum. It is interesting that 
also Potll, a dual with spectrum of type m ~ n, also leads to the simple thermodynamics 
in Fig. 8. 

4.7 The conformal window 

A detailed picture of thermodynamics in the conformal window is shown in Fig. 21. Here 
p/T 4 , i.e. the effective number of degrees of freedom, is plotted for some values of Xf > x c . 
At large T it is normalised so that it approaches the SB limit (3.38) for any xj. For T 
approaching zero, p/T 4 approaches another constant, the value of which decreases when 
Xf approaches x c from above. For all Xf, the vacuum phase has zero pressure, and at the 
limit T — > there is a transition from the black hole to the thermal gas phase. When Xf 
approaches the upper end of the conformal window 11/2, the behavior of the curves can 
be worked out analytically in perturbation theory [46] since the coupling then is small. 

For the present potential the finite temperature transition between the low and high 
temperature phases inside the conformal window is a smooth crossover. Fig. 21 also plots 
the interaction measure, the maximum of which defines the critical temperature for this 
crossover. Note that even if the transition here is smooth crossover, the transition can 
also be of 1st order in different theories [33]. What determines this behavior is the overall 
magnitude of the beta function. For illustration, consider the beta function of large Nf 
QCD, f3 = —boX 2 — b\X 3 + . . . . The values of the coefficients behave as bo ~ 0(1) -C \b±\, 
while the results of [33] suggest 6o ~ 0(10) ~ for 1st order phase transition. For the 
models we have considered here, we find that the nonperturbative beta function extracted 
from the gravity solution is small over the entire range < A < A* inside the conformal 
window. 
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Figure 21: Thermodynamics for some values of Xf within the conformal window, computed for 
Potl. Note that in the conformal window always r = and the functions a(A), n(X) do not affect 
the result. 



The large temperature values appearing in Fig. 21 may appear somewhat surprising. 
However, they have a simple explanation. The region in which p/T is nearly constant and 
approaching its large T limit is the perturbative region A small. The conformal window 
is within < A < A*(x c ) and the upper limit is always small, <C 1, so that the conformal 
window is perturbative, down to T = 0. From Fig. 2 one sees quantitatively how this holds 
even somewhat below the conformal window. To 1-loop logT = exp(l/(&oA)) and clearly 
for bo — > this grows fast. Somewhat more quantitatively, the beginning of the large 
T region can be estimated by computing the value of T for which the 2-loop correction 
term in the perturbative expansion of l/A(/i) equals the 1-loop term. One finds that 
the 2-loop correction is smaller than the 1-loop term if T > (2logT)^ bl ^ b o, |&i|/6g = 
3|13xj — 34| / (2(11 — 2xf) 2 ). This is always true for T > 1 = A if Xf < 3.6. However, for 
Xf > 3.6 this gives a lower limit of T which grows extremely fast when xf grows within the 
conformal window. T crossover is somewhat below the solution of this equation. Numerical 
values are in qualitative agreement with Fig. 21. 

4.8 The limits Xf — > and Xf — > x c 

The V-QCD models at Xf = are equivalent to an IHQCD model with potential V g . One 
thus expects that the hadronisation transition Th will approach the 1st order deconfining 
transition of SU(iV c ) YM theory 16 when xt — > 0. However, for Potl and large Wo this 
cannot be the case, since the hadronisation transition does not exist for very low xt, as we 
pointed out above. Indeed we see from the phase diagrams of Fig. 18, that the symmetric 
transition T s , which only exists for Xf ^ 0.4, is the precursor of the YM transition in this 
case. 

Let us then discuss in detail what happens in the xj — > limit in the two cases and 
at finite temperature. Since thermodynamics is determined by the set of T(\h) curves in 
Fig. 7, one should study how this configuration behaves when Xf — > 0. The T(A/ 1 ,t = 0) 
curve (shown in red in Fig. 7) exists only for A < A* and, according to the definition (2.40) 

16 Note that strictly speaking the limit of YM theory demands Nf = and falls outside the Veneziano 
limit Nf — > oo of QCD. This may explain the nontrivial structures observed at Xf — > 0. 
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Figure 22: T(Xu) for various small values of Xf and for potential I, Wo = 24/11 (Left) or for 
potential II, Wo = 12/11 (Right). The black curve is the IHQCD limit. The chirally unbroken 
T u (Xh) = T(Xh,r = 0) branch asymptotes to the IHQCD curve as Xf — >■ 0, for both potentials. 
The chirally broken Ti,(Xh) = T(Xh,Th(Xh,m q = 0)) branches behave very differently for PotI and 
Potll. For PotI T\, is absent at such low Xf and all phases are chirally symmetric (see also Fig. 18). 
For Potll the curves follow very closely T s and, correspondingly, « T s (see Fig. 15). 



A* — > oo when xt — > 0. In more detail, the limit is given by 



A * /3 = Wp V s V(K) = ^^A* /3 7mX. (4.3) 

where the equation on the left determines A* while the equality on the right expresses 
V(A*) as a function of A*. Eq. (4.3) is valid both for potentials I and II. The red curves of 
Fig. 7, therefore, stretch to the right when Xf — > 0. Quantitatively what happens is shown 
in Fig. 22 and one sees that they approach the T(Xf l ) curve of IHQCD when Xf — > 0. This 
is as expected since only V g remains in the limit. It is thus obvious that T s approaches 
the transition temperature of IHQCD as Xf — > (but it may be a transition between two 
metastable phases rather than a physical transition). 

To find the relative behavior of T s and one needs the asymptotic properties of 
the curves Tb(Xh) = T(Xh,Th(Xh,rriq = 0)) (shown in blue in Fig. 7) which only exist for 
A > A cn d. Here PotI and Potll behave in considerably different ways, as is already seen 
from Figs. 15 and 18. 

The crucial difference between PotI (at large Wo) and Potll comes from the fact 
that for PotI the value of A en d (the endpoint of the blue Tb(Xh) curves in Fig. 22) grows 
rapidly when Xf decreases, while for Potll A en d remains almost constant. Since always 
T u (X cn d) = Tb(X cn d) and T u (Xh) decreases rapidly at large A^, also the temperature T en d 
becomes small at small Xf for PotI. This drives the whole curve T^A/J towards zero and 
since is determined by integration along Tb(Xh) also — > 0. Finally A cnc j ceases to exist 
when Xf goes below a critical value ~ 0.25, the temperatures T en( j and reach zero, and 
the low temperature chiral symmetry breaking phase disappears. 

For Potll Tj,(Xh) follows very closely Tthqcd above it (Fig. 22) and it is thus natural 
that T/j ^ T s and that they approach the same limit. T s is actually metastable (Fig. 15). 
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Figure 23: The temperatures T u (Xh), Tb(Xh) of the two black hole branches, unbroken and broken, 
scaled to the temperature T ent j for Potll, W SB normalized. The values of Xf from top to bottom 
are 2.5, 2.75, 3, 3.3, and 3.5. Compare to Figs. 9 and 13. 



One can also illustrate the connection of the behavior of A en d to the BF bound of 
the tachyon (Fig. 3) by analyzing the linearized tachyon equation motion as discussed in 
Appendix D. 

In the limit of Xf — > x c one expects that all dimensionful quantities sensitive to the IR 
vanish as specified by Miransky scaling (1.8). All our numerical results are compatible with 
this, but conclusive numerical verification would require extensive further work. Analytic 
arguments supporting the scaling, similar to those presented in Sec. 10 of [22], can also be 
constructed in the finite temperature case. We shall here, however, only briefly comment 
on the scaling as well as the overall behavior of the solutions 0jS Xf r Xq from below. 

We start with the case of zero temperature which is simpler. For Xf < x c the dominant 
background is the one with nontrivial tachyon, and chiral symmetry is broken. As Xf — > x c 
the solution comes closer and closer to the fixed point, and the near conformal region grows. 
As it turns out, the pieces for A > A* and A < A* approach separately fixed solutions in 
this limit, which do not talk to each other. 17 Thus any observable which can be expressed 
only in terms of either the UV or the IR solution approaches a fixed value in the x f — > x c 
limit. The ratio of the characteristic scales of the two pieces diverge as specified by the 
Miransky scaling factor of Eq. (1.8). 

It is hard to find simple examples of such observables at zero temperature, but as 
it turns out, at finite temperature there are plenty. However, the analysis of the limit 
is more involved, since we have the additional parameter A^ we can be either in the UV 
(< A*) or in the IR (> A*) regions. The curve which determines the main features of the 
thermodynamics is Tb(Xh), which lies mostly in the IR region. Its endpoint A en d is however 
smaller than A*. When A^ > A* we expect that the background solution breaks into two 
parts similarly as for T = in the limit Xf — > x c , and the temperature is determined solely 

17 More precisely, keeping fixed the scale Ajj^ defined by the IR expansions, the background approaches 
pointwise a "IR" limiting solution which flows from the good IR singularity to the fixed point at A = A*. 
This solution which approaches A, from the "wrong" side is possible due to the presence of the tachyon. 
Keeping the UV scale Ayy fixed instead, the "UV" limiting solution is the one that flows from the IR fixed 
point at A = A* to the standard UV fixed point at A = with zero tachyon. 
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Figure 24: Examples of equation of state with qq condensate. Left: Type II SB-normalised 
potentials at Xf = 3 (compare Figs. 9 and 13). Note the scaling of a by a factor 100. Right: Type 
I potentials with Wo = 12/11 and Xf = 3.3 (compare Fig. 18, middle panel). 



by the IR piece. In the limit xj — > x c we find that A cn d — > A* from below. Therefore 
the whole Tj,(Xh) curve is in the IR region in the strict Xf — > x c limit, and it is plausible 
that it takes a fixed shape. 

This behavior is supported by the numerical study of Fig. 23, where we plot the 
temperatures of the two black hole branches as functions of Xh for Potll with Wq SB 
normalized. The xy-dependence of the curve T^A/J is, up to the overall normalization, so 
small for Xh 3> A* that it cannot be resolved from the plot even at relatively low xt < 3. 
The main effect with increasing Xf is that the Tb(Xh) curve is visible down to lower and 
lower Xh as A* decreases slowly, while the shape of the curve remains fixed. The curve 
T u {Xh) approaches a vertical line when scaled to T en d, reflecting the fact that A en d — > A*. 

The values of all the critical temperatures (except T crossover ), as well as all thermody- 
namics up to the transitions, are determined by T^A/J as Xf — > x c . Therefore, we expect 
that the thermodynamics "freezes" in this limit, in the sense that all ratios of the critical 
temperatures approach fixed values. Moreover, the parts of the thermodynamical func- 
tions which are determined by the IR solutions, are expected to have well defined limits. 
While we have not proven these statements, they are strongly supported by the numerical 
study of Fig. 23. Notice however, that our data only extends up to xj = 3.5 which is still 
well below the critical value x c ~ 3.70. Therefore we cannot exclude the possibility that 
something drastic happens for Xf even closer to x c . 



There is a subtlety here as comparing temperatures at different Ah requires fixing the units of energy, 
which we do by calculating Atjv- However for Ah > A* the UV part of the solution, and hence practically 
the units of energy, become independent of Ah in the limit xf — > x c . Thus we could equally well define the 
units in terms of the behavior of the IR piece of the solution near A = A, . 

19 This is observed numerically, and can be understood by studying the violation of the BF bound in the 
spirit of Appendix D. 
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Figure 25: The condensate for the Pot I, Wq = II phase diagram in the right hand column of 
Fig. 18 at Xf = 0.3 in linear (left) and log (right) scale. Note the scaling of a by a factor 1/250 in 
the left panel. The right panel shows all the three transitions in Fig. 18. 



5. The chiral condensate 

In principle, it should be straightforward to extract the chiral condensate oc a from the 
tachyon solution in the UV 

r(r) ~ or (-ln(Ar)) 3 /( 2fc(, \ r^0, (5.1) 

as the quark mass is set to zero. However, in this model the task is actually very demanding 
due to the logarithmic corrections (i.e., the running of the condensate) and the fact that 
the numerical solutions have a tiny residual quark mass due to limited numerical precision. 
These issues and their resolution are discussed in Appendix E. 

Examples of the correlation of a with that of the free energy are shown Fig. 24. One 
expects that a jumps in a 1st order transition. The first case is thermodynamics for SB- 
normalised type II potentials at Xf = 3, also studied in Figs. 9 and 13. Cooling from large 
T chiral symmetry first breaks at the 2nd order transition T en( j and the condensate starts 
from zero and increases with further cooling. At Th the system experiences a 1st order 
transition and a jumps the amount shown in the figure. Below that a remains constant in 
the present models, which does not describe the thermodynamics of the low T phase. In 
that the degrees of freedom are Nj massless Goldstone bosons. 

A second example is thermodynamics with condensate of Type I potentials with Wq = 
12/11 and Xf = 3.3. This case is special in that in it the 1st order line splits in two 1st 
order transitions at T/j, T\2 if Xf ^ 2.8, as shown explicitly in Fig. 18, middle panel. Again 
the highest temperature transition is a 2nd order one at T en( }, at which the condensate starts 
to grow when the system is cooled. The condensate grows up to the value a = 4.537A 3 . 
Then there is a 1st order transition at T\2 with a jump in a and latent heat: 

Act Ae , , 

= 0.508, -3- = 2.29 (5.2) 
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Figure 26: The condensate as a function of T en( ± — T in log-log scale. The dots are our data, 
extracted from the background solutions, and the lines are fits having the form a = C\fT eni ^ — T. 
The choices of potentials are, from top to bottom: PotI with Xf = 3.5 and Wq = 12/11 (red); PotI 
with Xf = 0.3 and Wq = 24/11 (blue); and Potll with xt = 3 and SB normalized Wq (magenta). 





and finally a very weak transition at with the value 

~ = 1.03-1(T 7 . (5.3) 

It is clear that the value of a also jumps at the latter transition, but the size of the jump 
is so small that we could not extract it reliably from the numerics. 

As a third example, consider the case Pot I, Wq = at Xf = 0.3, which is very special 
in that is very small and it is T s which dominates, as is seen in the right hand column 
of Fig. 18. The magnitudes vary so much that all the transitions can be presented only on 
log scale, see Fig. 25. 

At the 2nd order transitions the condensate goes to zero continuously as the temper- 
ature approaches T en d from below. The curves in Figs. 24 and 25 seem to be compatible 
with the standard expectation a ~ V^cnd — T. We study this more precisely in Fig. 26 
where we plot our data for the condensate for T close to T en d in the log-log scale for various 
choices of the potentials. The data are compared to the lines a = C\/T cn d — T, with the 
coefficients C chosen such that the lines overlap with the tails of the data at small r cnc j — T. 
It is convincing that 1/2 is indeed the correct exponent. 

6. Nonzero quark mass and thermodynamics 

We have so far mainly discussed the case of vanishing quark mass, for which chiral symmetry 
may hold. Effects of nonzero quark mass have been mentioned in subsection 3.3.2 and 
described in Figs. 6 and 7. They follow from the fact that tachyonfree black hole solutions 
do not exist, as a reflection of the disappearance of chiral symmetry. However, numerically 
the tachyonic small-m 9 solutions are very close to the zero mass tachyonfree solutions in 
the UV at small A, as is seen by comparing left and right panels of Fig. 7. Thus chiral 
symmetry is always broken even in the high-T phase, but quantitative effects are small at 
large T. 
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The effects of small nonzero m q are shown quantitatively in Fig. 27 and can be sum- 
marised as follows: 

• The main effect follows from the fact that for nonzero m„ even the high T phase is 
chirally broken. The curve marked m q = is the same r^(xj) as that in Fig. 13. For 
this case the phase at T > Th is chirally symmetric and chiral symmetry is broken 
when T decreases below T^. For nonzero m q also the phase at T > Th is chirally 
broken and the effective order parameter of the transition is the jump in entropy or 
energy density. There is also a jump in the condensate, but the condensate is nonzero 
also for T > Th{xf). 

• For Xf clearly below x c the effects of small m q on the phase diagram are small. 
Particularly interesting is the pattern of approach towards m q = 0. The smaller m q , 
the higher is the value of Xf where the curves start deviating significantly. 

• The 2nd order transition at T enc [ becomes a continuous one. This is obvious from 
Fig. 7, there is no discontinuity near T en( j. However, a remnant of the genuine tran- 
sition is a maximum of interaction measure, also plotted in Fig. 27. 

• At large Xf the conformal window and Miransky scaling disappear. For m q = the 
transition temperature vanishes when x c — Xf — > as dictated by Miransky scaling. 
The smallest nonzero mass destroys this effect and T/ 4 curves upwards towards larger 
values. 

• The effect of nonzero mass could also be seen by plotting the beta function for values 
of Xf within the conformal window. For m q = only the tachyonless solutions matter 
and they extend only up to A* in Fig. 7. The beta function /3(A) only exists for A < A* 
and (3(X*) = at the IR fixed point A*. For m q > the beta function comes close to 
A* but continues past it to larger values of A. 

7. Outlook 

In this paper, we have used bottom-up holography to study the thermodynamics of models 
that are in the universality class of QCD with massless quarks in the Veneziano limit, (large 
N c and N f but fixed x f = N f /N c ). 

The temperature dependence of the pressure p(T) and the condensate cr(T) = (qq)(T) 
was computed at various x/ up to the loss of asymptotic freedom at Xf = 11/2 with 
a conformal window appearing at x c < Xf < 11/2, x c ~ 4. The singularities of these 
quantities define the phase diagram of the system. 

One expects that the system has two phases, a low temperature phase with broken 
chiral symmetry and a chirally symmetric high T phase. The simplest alternative is that 
these are separated by a 1st order transition at some T = T(xf), see Fig. 8. In holography, 
the high T phase is a metric with a black hole and a new feature appears: this phase can 
be either chirally symmetric (no tachyon) or chirally broken (nonzero tachyon). The same 
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Figure 27: The behavior of Th(xf) (in units of l/Cyy, over the Xf range in the figure this unit 
changes by about 30%) for m q — and for small values of m q . For m q = the dashed line shows 
the true 2nd order chiral symmetry breaking transition. For nonzero m q the line shows the position 
of one maximum of the interaction measure, a second one gives the usual T crossov er (which is almost 
independent of m q ). 

doubling applies to the low T thermal gas phase. The phase structure can be correspond- 
ingly more complicated. For given gluonic and fermionic potentials of the V-QCD action 
the thermodynamics is fixed and calculable. 

A typical prediction of the model is that there indeed is a 1st order line for Xf from 
up to some value < x c . In holography this is a transition at some T = T^^Xf) between a 
non-tachyonic black hole metric and a thermal gas metric with a tachyon. In field theory 
language one would say that at this temperature there is deconfinement and chiral symme- 
try breaking with coinciding deconfinement and chiral temperatures, Td = T x . However, 
a new interesting feature is that at larger Xf this 1st order line can split in two: first at 
higher T chiral symmetry is broken in a 2nd order transition, then at a lower T there is a 
1st order deconfinement transition (see, e.g., Fig. 13), T x > Td- 

A particular feature of the phase diagram is that with x f approaching x c all transition 
temperatures, as all mass scales, decrease as specified by Miransky scaling. Associated 
with this approach is quasiconformality and walking. At very large T there is always a 
weak coupling region into which one enters at T crossover . This increases when Xf grows, 
behaves regularly at x c and continues further into the conformal window increasing faster 
and faster (see again Fig. 13). 

Detailed predictions of the model depend on the gluonic and fermionic scalar poten- 
tials. There are several physical constraints in deriving them, but they are not uniquely 
determined. One crucial constraint is missing: in contrast to the application of this model 
to hot SU(A r c ) theory [27], there is no 4d lattice data in the large N c , Nf limit. Thus 
one is genuinely making predictions and we find that also phase diagrams deviating rather 
drastically from the above baseline prediction follow from holography. Particularly striking 
examples are shown in Fig. 18. In these, the chirally symmetric phase can extend all the 
way down to zero temperature. 

Although we do not expect this to happen in QCD, it may happen in other large N- 
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theories. The last diagram of Fig. 18 has a structure that is very reminiscent of the phase 
diagram of high-T c cuprates, although here we are at zero charge density. In particular 
the intermediate dome-like structure with chiral symmetry breaking corresponds to the 
superconducting dome [50] in the cuprate diagram. 

One may also note that the models predict a very rich structure of metastable states. 

There are many directions in which the model could be further sharpened and devel- 
oped: 

• All the computations are numerically very demanding and much further work would 
be useful. This holds, in particular, for the approach to the conformal window, 
Xf — > x c , and for computations of the condensate. 

• Effects of non-zero quark mass on the thermodynamics should be studied beyond the 
discussion in Section 6. 

• The extension to Xfi flavors of mass mi, etc., is possible but requires a non-abelian 
version of the Sen action. This is in principle possible to deal with, although we have 
much less information from string theory on the details of such an action. 

• It would be very interesting to accumulate lattice Monte Carlo data in the Veneziano 
limit. For pure SU(iV c ) theory already N c = 3 is very close to N c = oo [51, 48]. 
One might thus see that, e.g., Nf = 4,8,12 at N c = 3 would already give useful 
information. 

• A chemical potential for baryons should be included. This necessitates the inclusion 
of a bulk U(l) baryon vector field A^. It would be particularly interesting to know 
the fate of the T x > phase at large Xf and fi. 
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APPENDIX 

A. UV expansions 

The expansions near the standard UV boundary can be computed in a straightforward 
manner. Let us fix the location of the boundary at r = 0. The blackening factor / is very 
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close to unity in the UV, 



f(r) = l + 0{r 4 



(A.l) 



Therefore the leading finite temperature and zero temperature expansions of the various 
fields coincide. Moreover, as the tachyon vanishes at least linearly in r for r — > 0, it can 
be set to zero when solving for the leading behavior of the coupling A and the warp factor 
A. Hence the expansions take the form familiar from earlier works [26, 27]. We reproduce 
here the leading expansions of A and A as well as the expansions of the non-normalizable 
and normalizable tachyon expansions both in terms of r and A. 

A.l Fields A and A 

Setting the tachyon to zero, the equations of motion for A and A involve the effective 
potential 



V cS (\) = V g (\)-x f Vf(\,0) 



12 



[1 + VxX + V 2 X 2 + 



(A.2) 



Then the (leading) UV expansions of A and A can be found by substituting suitable Ansatze 
in the equations of motion. The result reads 



A(r) 
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£uv 
95 



+ 



91n(rA) 
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+ iln[-ln(rA)] 



-23 + 



64V 2 



In 



ViA(r) 



ln(rA) 2 
ln(rA)] 



+ o 
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ln(rA) ; 



(A.4) 



Two combinations of the series coefficients of the effective potential appear here. As the 
potential is matched with perturbative QCD, they become 
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(A.6) 



where bi are the coefficients of the perturbative QCD beta function. Notice that these 
coefficients are indeed the same for all potentials used in our study and in particular 
independent of the choice of Wq. 

Let us also present the expansions in terms of A, as we use it as a coordinate in all 
numerical calculations. The result after the conversion reads 

4 

9A 



lnr(A) = -^ + ln(£ uv ) 



(A.7) 
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A. 2 The tachyon 

As the tachyon is decoupled near the UV boundary, its UV behavior can be studied by 
inserting the expansions calculated above for A and A into the tachyon EoM. We also 
develop the potentials as series in the UV: 

12 



Kff(A) = V g (X) - x f Vf(X, 0) = -f- [1 + ViA + V 2 X 2 + 

^uv 



xVf(X) = W + WtX + W 2 X 2 + 



k{X) _ 2C 



uv 



[1 + kiA + k 2 A 2 + 



(A.9) 



(A.10) 



Here the leading coefficient of n/a was already fixed in order to have the correct UV mass 
of the tachyon [29]. It is enough to study the linear terms in the tachyon EoM, which 
become 
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Here matching with the perturbative anomalous dimension of the quark mass in QCD gives 
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where 70 is the leading coefficient of the anomalous dimension of the quark mass in QCD. 
The result can be again written in terms of A, and it becomes 
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+ aCfj V e- 3A (lnA)^ [l + O (A- 1 )] 

A. 3 Finite temperature 

The basic relations 
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can be evaluated in the UV by inserting from (A. 3) 
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Terms of the order of In 2 (In) /In 3 are neglected; for these, see [52]. One finds 
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For the quantity bh/T needed for the latent heat one has 
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B. IR expansions 



B.l Zero temperature 

Here we first discuss the expansions near the "good" IR singularity at zero temperature. It 
is the particular solution which can be lifted to finite temperature. In the IR, the tachyon 
potential in the DBI action is expected to be exponentially suppressed. Therefore the 
tachyon is again decoupled, and the IR behavior of A and A can be solved separately from 
that of the tachyon. Moreover, the IR expansions of A and A are exactly the same as in 
IHQCD. We will anyhow repeat the discussion for the particular asymptotics of V g that 
matches well with the IR properties of QCD [26], which covers all potentials in this article. 

B.l.l A and A 

Let us assume that the potential V g has the asymptotic behavior 

V\ v 2 
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Then the asymptotic solution reads 
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20 Ref. [26], second paper, equation (D.3), has a different constant in the expansion of T. We have checked 
the 1/3 here also numerically. 
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The IR scale R = 1/Afr is an integration constant here. 21 Recall that r does not appear 
explicitly in the equations of motion, and therefore there is also an integration constant 
related to r: we have the freedom of shifting any solution by r — > r + Sr. The solution 
having the simple r-dependence of Eqs. (B.2) and (B.3) corresponds to a special choice of 
Sr. It will have its UV boundary at an arbitrary value of r (rather than at r = 0). If Sr is 
fixed instead by requiring the UV boundary to lie at r = 0, a corresponding shift must be 
added to the asymptotic formulas. For our choice of V g , 
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if we set Co = 1. 

Using A as the coordinate, the result reads 
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Various other combinations may be useful. In thermodynamics one needs b = e A in 



terms of A; from (B.9) one can invert: 
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For the pair of functions q(A) = e r'(A), X(A) used in numerics one can derive, for A 
— oo, 
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B.1.2 The tachyon 

The IR expansion of the tachyon depends on the large- A asymptotics of the potentials Vf, 
a, and k. Recall that the tachyon potential V/(A,r) needs to vanish in the IR [29] in order 
to have correct kind of flavor anomalies. All power-law asymptotics for the potentials were 
analyzed in [22], and two different acceptable cases were chosen as examples. They are: 



I Asymptotics with 



a(A) ~ A ; k(A) ~ A 4 / 3 



Vfo(X) ~ A T 



(B.13) 



21 It is not independent of the scale A = Ayy of the UV expansions, i.e., the complete solution from UV 
to the IR will have fixed Atjv/Ajj^. 
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where f < 10/3. This case includes the potentials I and I* of this article (for which 
f = 2). The tachyon diverges exponentially for r — > oo (A — > — oo), 



r ~ e C/ i ~ e Cl ^ (B.14) 
where the coefficient reads for potentials I 

= 81 3^115 -Ifa^ll-x,) 

812944 2V6 v ; 

II Asymptotics with 

a(A) ~ A* ; k(A) ~ A 4/3 ; V^A) ~ A f (B.16) 

where a > and f can take any value. This case includes the potentials II and II* 
of this article (for which a = 2/3 and f = 2). The tachyon diverges for r — > oo 
(A — > — oo) as 

r-C^^-C/K-^) 1 / 4 (B.17) 

where the coefficient reads for potentials II 
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B.2 Finite temperature 

We will work out the finite temperature IR expansions in A-coordinates. Instead of writing 
down the explicit expansions as above, it is more convenient to state the relations between 
the coefficients of the series expansions. We start by defining the series 



(B.19) 



f = ef' h + 0(e 2 ), f'(0) = f h + ef^ + O(e 2 ), 

q = q h + eq' h + 0(e 2 ), (B.20) 

A = X h + e \' h + 0(e 2 ), (B.21) 

T = Th + eT i h+ i e \ii + 0(e 3 ), (B.22) 

where e = A — A^ is the distance from the horizon, which lies at A = A^, and all coefficients 
are to be evaluated at the horizon. The key input here is f(Ah) = 0. Inserting to the 
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equations of motion one can solve for six of the nine coefficients listed above: 



1h 

ft 
% 
A', 



-4A + 4 
Jh 



(fD 2 



qh[i + f h h 



J M(dxV g -d x V f ), 

£1 ' 

fh K h 

9d T V f {A + B + C) + D 

i2K\vj{v f -v g Y ■ 



(B.23) 
(B.24) 
(B.25) 
(B.26) 
(B.27) 
(B.28) 



with the abbreviations 



A 
B 
C 
D 



6\ 2 h K' h v?(d x v a -d x v f) , 



v?[8d?v f - 3\ 2 h (d x v g - d x v f )(K h (d x v g 



3d x V f ) + 2 K ' h V g )}, 



-2V f [Qd T V f + V g (4%V f + 3X 2 h K h d x V f (8 x V g - d x v f ))}, 
27X 2 h K h d T d x V f Vf(V g - V f )(dxV g - d x Vf). 



(B.29) 



Here V g = V g (Xh), Vf = Vf(Xh,Th), n' h = dn(Xh)/dX, Kh = n(Xh)- The so far unspecified 
three coefficients Xh, Th, and ft remain as free parameters. However f' h will be fixed by 
requiring the standard normalization of the blackening factor / — > 1 in the UV. Therefore 
the physically relevant parameters are and Th, which can be mapped to the temperature 
and the quark mass after the full solution has been found. 



C. The quark mass and the Efimov solutions 

As detailed in [22], the existence of the Efimov vacua is tightly linked to the tachyon mass 
at the IR fixed point, plotted in Fig. 3. In particular, the existence of the full Efimov 
tower of vacua with arbitrary number of tachyon nodes is guaranteed if the tachyon mass 
violates the BF bound. The same holds at finite temperature: one can always tune Xh and 
Th such that the solution comes arbitrarily close to the fixed point. When the BF bound is 
violated, the tachyon solution is oscillatory in the vicinity of the fixed point. Thus, when 
approaching the fixed point the tachyon will achieve arbitrary many nodes, which signals 
the presence of the full Efimov tower. In this case the dependence of the quark mass on 
Xh and Th is the "standard" one, i.e., qualitatively as in Fig. 5. 

There are, however, some cases where either the fixed point is absent, which is the 
case for potentials I* and II* at low xf, or the BF bound is not violated, which is the case, 
interestingly, for potentials I at very low xj (as well as in the conformal window for all 
potentials). In such cases the picture can be different from Fig. 5. We shall not give a 
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detailed description of all possible cases here, but rather discuss some of the main features 
and give examples. 

The curve The (which actually starts at A*) exists if and only if there is a fixed point. 
If there is no fixed point, the solutions are expected to reach the standard UV boundary 
for all values of and r^. For the curves r^o, Thi, ■ ■ ■ the situation is more complicated. 
At least few of these curves may still exist even if there is no fixed point or if the BF 
bound is satisfied at the fixed point. Their existence at asymptotically large A^ is linked 
to the existence of Efimov solutions at zero temperature: taking A^ — > oo with fixed 
along the curves, the finite temperature Efimov configurations converge towards their zero 
temperature counterparts. In particular, we expect that the chiral symmetry is broken 
at zero temperature if and only if Tho exists at asymptotically large A^. We have found 
numerically that the curves are always absent in the conformal window, xf > x c , so that 
chiral symmetry is intact. This turns out to be the case also for potentials I at large Wo 
and low Xf, but only in a part of the region where BF bound is satisfied at the fixed point. 
See also the phase diagrams in Fig. 18 of Sec. 4 which show that chiral symmetry is intact 
at low Xf. For potentials I* and at low Xf, where no fixed point exists, the curves are 
also absent, and chiral symmetry is unbroken. In this case the m q = thermodynamics is 
determined by the r = solution and is qualitatively similar to the Yang-Mills one (see 
also Fig. 19). For potentials II* however, at least the leading solution Tho can always be 
found and chiral symmetry is thus broken at low temperatures (see Fig. 17). 



D. Computation of A en d 

One can also illustrate the connection of the behavior of A en d to the BF bound of the 
tachyon (Fig. 3), assuming that we have chosen a set of potentials and value of Xf such 
that the IR fixed point exists. First we recall that A em } can be defined as the endpoint 
of the t/jo(A^) curve which gives the (non-node) solution with nontrivial tachyon and zero 
quark mass (Fig. 5). In particular, as A^ approaches A en d from above, tends to zero, 
and we expect that the whole tachyon solution from the boundary to the horizon becomes 
small, and the tachyon decouples from the other fields. Therefore, in order to define A en d 
it is enough to study the behavior of the tachyon based on the linearized tachyon EoM, 
evaluated on a fixed background, obtained by setting the tachyon to zero. 
The linearized tachyon equation has the form 



where 



r"(r) + Fir'(r) + F 2 r(r) = (D.l) 



Ft = M'(r) + + A '(r)^W + X^rf^M^, (D.2) 
f(r)K(X) ■ 

Here A(r), A(r), and f(r) are the solutions of the EoMs for r = 0, which are the same for 
potentials I and II. The drastic difference between the potentials, as suggested by Fig. 3, 
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thus arises only through the appearances of a and k in the coefficients (D.2). The regular 
tachyon solution, which is finite in the IR, obeys 

^ = - Hm § (D.3) 

since the double-derivative term in (D.l) is negligible near the horizon. 

Nodes of the regular solution to the linear tachyon equation can then be used to 
determine A cn d. For small A/j perturbative analysis applies and it is not difficult to see 
that the solution is monotonia, without nodes. When A^ increases the equation becomes 
nontrivial and has to be studied numerically. Usually we observe, that beyond a critical 
value of Xh a tachyon node appears in the UV. The leading tachyon behavior in the UV is 
controlled by the quark mass, which has to vanish at the critical value. We thus identify 
the critical value as A en d> which was defined as the endpoint of the curve where m q = 0. 
Thus the regular solution to the linearized EoM has no nodes for A^ < A cn d and one or 
more nodes for A^ > A en d- It can also happen that A en d does not exist, and the tachyon 
nodes are absent for all A^. 

Since A^ can take values from zero to A*, we can construct backgrounds which get 
arbitrarily close to the IR fixed point at A = A*. If the BF bound for the tachyon is 
violated at the fixed point, the tachyon must have nodes as Xh — > A*. We can conclude 
that Aendi and thus also the curve Tho, exist in this case. This makes sense, since when the 
BF bound is violated, chiral symmetry breaking takes place also at zero temperature, which 
means the the curve tho exist also at asymptotically large A^ as discussed in Appendix C. 

We can also say something about A en d in the probe limit Xf — > 0. For Potll it seems 
that it approaches a fixed value as seen from Fig. 22 (right). This value can be found 
by solving the linearized tachyon EoM with a background evaluated at Xf = (i.e., the 
IHQCD solution), and by checking if a special value of A^ (identified as A en d) can be found 
where nodes emerge in the tachyon solution. Notice that A* goes to oo in the probe limit so 
that Xh can take any value. Existence of the limiting value of A cn d as xj — > thus requires 
that the tachyon has nodes in the limit Xh — > oo after first taking the probe limit. Since 
zero temperature solutions are obtained for Xh — > oo, it is plausible that the behavior of 
Fig. 22 (right) is seen if and only if the probe limit system admits tachyon solutions with 
nodes (in other words, chiral symmetry breaking) at zero temperature. Recall that for PotI, 
for which the different behavior of Fig. 22 (left) is found, chiral symmetry is unbroken at 

low Xf. 

E. Computation of the condensate 

In principle, the condensate for an m q = system could be computed from the UV expan- 
sion 

r ( r )/£ uv = (Tr 3 (-ln(Ar)) 3 /( 2bo \ r -> 0, (E.l) 



with 

A - In (A£ uv ) = r^r + % ln(M(A)) = - In (Ar) , A -> oo, 



(E.2) 
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where we dropped corrections of 0(A ). Using this one can define 

hxa(A) = lnr(A) - ln£ uv + t-^— + ^MM(^)) + Hb X(A)), (E.3) 

which approaches In a for A — > oo. 

However, our solution for the tachyon, which is obtained numerically by shooting from 
the IR, will have a linear term r ~ m q r with a tiny quark mass (typically m q ~ 10 -7 ), 
because the IR boundary conditions cannot be fine tuned beyond the numerical accuracy 
of the code. The linear term will dominate over the cubic one of Eq. (E.2) in the deep 
UV. In order to calculate the condensate, we need to separate the linear and cubic terms 
from the numerically computed t{A), and use the cubic solution in Eq. (E.3). For cr(A) 
to be a good approximation to the condensate a, we need to have A ~ hundreds. Direct 
separation of the linear m q term in this region requires numerical accuracy on the level of 
e -hundrcds^ -^^j^ j s practically impossible to achieve. 

To illustrate the difficulty and its resolution, consider a concrete case. Let us take 
Potential II, SB normalised, £tjv = (1 + \ x f) 1 ^-, %f = 3. This system, when cooled, 
has a 2nd order transition at T en( j = 1.1587X, above a 1st order transition at 7\. This is 
concretely seen in Fig. 9. Since chiral symmetry is broken at T en d we expect that c(T) 
starts growing from zero at T en d and grows when the system is cooled towards T^. As an 
example, we evaluate the condensate when T has been cooled to T = 0.95T cn d = 1.17ft. 

Numerical solution of Einstein's equations required knowing the values of A^, leading 
to a certain T with m q = 0. For this potential and T they were A^ = 0.4017564, = 
Tho(Ah,m g = 0) = 0.217984. The computed r(A) is shown in Fig. 28. For A up to about 
10 one discerns the required r 3 ~ e~ 3A behavior, but beyond that r ~ e~ A sets in and 
extends up to the end point of the computation at A = 400. It is impossible to shoot from 
the horizon and get m q = more accurately; note that the tachyon has already decreased 
to 10~ 14 from 0.22 at the horizon. 

To impose r(r) ~ r 3 one must shoot from the boundary, r = 0, A = oo. In this limit 
the evolution of r decouples from the other bulk fields, of which only X(A) is relevant since 
/ ~ 1. We can thus integrate the tachyon equation from some large A(= 400) using the 
\(A) from Einstein's equations and imposing as the initial condition t(A) = e~ 3A , t'(A) = 
—3e~ iA with small enough normalisation. 22 The result is plotted as the curve ttjv(^4) 
in Fig. 28. One observes that in the range A = 2. ..10 the curve behaves accurately as 
a constant xtir (A) and the normalisation can thus be determined. In this way the true 
t(A; m q = 0) plotted in Fig. 28 is obtained. 

Now that the accurate t(A) is known, lna(A) can be plotted using Eq. (E.3), see 
Fig. 28. For the extrapolation it is even more convenient to plot as a function of X(A), 
see also Fig. 28. One obtains a nice linear behavior with the asymptotic value lnu = 
-5.1558, cr(T = 0.95T end ) = 0.005766. 



22 It is not important to have precisely correct UV boundary conditions, since corrections to the r ~ r 3 
solution will decay fast as the system is solved toward the IR. One should only make sure that the tachyon 
is much less than one in the whole region of interest (A > 0) in order to suppress nonlinear effects, or 
alternatively use explicitly linearized differential equation for the tachyon. 
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Figure 28: The tachyon computed for T = 0.95T cnc j = l.lT^. The curve i\ji{A,m q — tiny) is 
obtained by integrating Einstein's equations from the horizon and tuning m q = as accurately as 
possible. The curve r\jy(A) is obtained by integrating the tachyon equation of motion from the UV 
at A = 400 using the bulk field X(A) from Einstein's equation and imposing r ~ r 3 in the UV. The 
normalisation can be fixed by matching to T[r in the A = 2. ..10 range and a reliable t(A, m q = 0) 
for the true m q = tachyon is obtained. 




Figure 29: Extrapolation of the value of lncr to r = 0, A = oo, A = 0, for the potential and 
temperature as in Fig. 28 

If one used the original t(A, m q = 0) at the largest value of A, A = 10, where the 
r 3 behavior was obtained, one would have a = 0.0106. This is too large by a factor 1.82, 
not very far off, but actually slightly larger than the expected 10% error from neglecting 
the 0(A~ 1 ) corrections in Eqs. (E.2) and (E.3) at this value of A. If we tried using the 
t(A, m q = 0) solution directly, reliable extraction of a would thus require much higher 
numerical precision, as already mentioned above. 

After applying the procedure discussed above, the dominant error in the value of a 
arises actually from the matching of the two tachyon solutions that were obtained by 
shooting from the UV and from the IR. The solutions are not exactly proportional for 
< A < 10 due to nonlinearities in the tachyon EoM and coupling to other fields. The 
error can be reduced by introducing a further subtraction trick that effectively reduces the 
value of m q of the solution that was obtained by shooting form the IR, so that the matching 
can be done for slightly higher values of A where the coupling effects are considerably 
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reduced. 



We follow [22] and construct two solutions n t 2 with small but different values m q \ iq 2. 
Optimal choice is to take \m q \ as small as possible and choose one solution with a positive 
value and another with a negative one. Then we construct 



where the ratio m q i/m q 2 can be accurately determined as the ratio of the solutions t\/t2 
at large A (say A = 400). 23 The point is that the constructed tir, has its m„ several 
orders of magnitude closer to zero than either of the solutions n,2- Moreover, the residual 
dependence of a on m q is drastically reduced: the linear corrections cancel in (E.4) (see 
[22]). The improved value of a can now be found by matching tir with the solution ruv, 
which was obtained by shooting from the UV, as discussed above. 
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